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Semi-inclusive tt"*" electroproduction on protons has been measured with the CLAS detector at 
JefFcrson Lab. The measurement was performed on a liquid-hydrogen target using a 5.75 GoV elec- 
tron beam. The complete five-fold difi'erential cross sections were measured over a wide kinematic 
range including the complete range of azimuthal angles between hadronic and leptonic planes, 0, 
enabling us to separate the ^-dependent terms. Our measurements of <?!)-independent term of the 
cross section at low Bjorken x were found to be in fairly good agreement with pQCD calculations. 
Indeed, the conventional current fragmentation calculation can account for almost all of the ob- 
served cross section, even at small tt"*" momentum. The measured center-of-momentum spectra are 
in qualitative agreement with high energy data, which suggests a surprising numerical similarity 
between the spectator diquark fragmentation in the present reaction and the anti-quark fragmen- 
tation measured in e^e~ collisions. We have observed that the two ^-dependent terms of the cross 
section are small. Within our precision the cos2(/i term is compatible with zero, except for low-z 
region, and the measured cos term is much smaller in magnitude than the sum of the Cahn and 
Berger effects. 

PACS numbers: 12.38.Cy, 12.38.Lg, 12.38.Qk, 13.60.Hb 



I. INTRODUCTION 



The semi-inclusive leptoproduction of hadrons off the nucleon, eN — > e'hX, is an important tool allowing to 

study simultaneously the internal structure of the target nucleon and hadron creation mechanism. In Deep Inelastic 
Scattering (DIS) regime the semi-inclusive leptoproduction of hadrons can be described by perturbative quantum 
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chromodynamics (pQCD) combining non-perturbative distribution/fragmentation functions. Semi-Inclusive lepto- 
production of hadrons in DIS (SIDIS) can occur through current or target fragmentation [Ij] (see Fig. [T|). Current 
fragmentation is the hadronization of the struck quark, while target fragmentation is hadronization of the specta- 
tor. Both non-perturbative, soft fragmentation mechanisms factorizc from the hard virtual-photon/parton scattering 
amplitude in pQCD (see Ref. [1] for the current fragmentation and Ref. [1,01 for the target fragmentation). 

Inclusive Icpton scattering off the nuclcon and hadron production in e^e~ collisions allow one to study separately the 
fractional momentum dependence of the parton distribution functions for the nucleon and the parton fragmentation 
functions, respectively. The leptoproduction of hadrons in the current fragmentation region combines these two 
and provides additional information about hadronization and nucleon structure. In fact, for DIS, semi-inclusive 
measurements provide new information about the Transverse Momentum Distribution (TMD) of partons, which 
is important for understanding the role of orbital angular momenta of quarks and gluons [sf. Furthermore, the 
detection of a hadron in SIDIS introduces a flavor selectivity for the observed parton distributions. In contrast, target 
fragmentation is described by fracture functions, present only in the semi-inclusive reactions. 

The finite transverse momentum of partons in the initial state leads to an azimuthal variation in the cross section, 
as does the transverse spin of partons in the unpolarized nucleon. 

In order to achieve SIDIS regime sufficiently high beam energy is mandatory. With decreasing the beam energy 
higher order (in pQCD) and higher twist effects appear, spoiling the agreement between the experimental data and 
theoretical pQCD expectations. Therefore, only a comparison between the actual data and theoretical calculations 
can reveal the dominance of SIDIS dynamics in the experiment. Though, a good agreement between data and theory 
in one observable does not necessarily guarantee SIDIS dominance in others. 




FIG. 1: Schematic representation of current (a) and target (b) fragmentation processes in the virtual photon-proton center of 
momentum frame neglecting transverse momenta of particles. Blobs represent non-perturbative functions: the parton density 
function fi{x), parton fragmentation function D'l{z) and fracture function A{i{x,z) as given in Eg. 1291 Xi and X2 indicate 
two components of the undetected final state hadronic system X. 

Previous measurements 0, 0, H, have verified these factorizations experimentally and have tested pQCD 

predictions. Measurements of unpolarized semi-inclusive lepton-nucleon scattering have been performed at several 
facilities such as CERN (EMC Q), Fermilab (E655 0]), DESY (HI [8], ZEUS [|, HERMES [l^), SLAC pi, 
Cornell [H, [H, [3 and Jefferson Lab (Hall C) [3, • The last two experiments covered a kinematical region similar 
to the present measurement. Despite all of these measurements, open questions remain about the target fragmentation 
mechanism and the physics behind the azimuthal distributions. The measurements at high beam energies (EMC, 
E655, HI and ZEUS) covered a broad kinematic range, but lacked particle identification and the statistics to look at 
differential cross sections in more than two kinematic variables (the latter applies also to HERMES). Experiments at 
lower energies (SLAC, Cornell and Hall C of Jefferson Lab) using classical spectrometers measured cross sections only 
at a few kinematic points. To improve the current knowledge of semi-inclusive lepton-nucleon scattering one has to 
combine the broad coverage of high energy experiments with high luminosity and particle identification in order to 
measure the fully differential cross section for a specified hadron. 

Semi-inclusive hadron electroproduction, 7^(5) -I- p{P) h{ph) -I- X, is completely described by a set of five 
kinematic variables. The variables g, P and pu in parentheses denote four- momenta of the virtual photon 7„, the 
proton p and the observed hadron h. The letter X denotes the unobserved particles in the reaction. In this article 
we have chosen a commonly used set of independent variables: the virtual photon four-momentum transfer squared 

= —q^ AEqE' sm^ |, the Bjorken scaling variable x — — 2Mi7' ^^"^ virtual photon energy fraction 

carried by the hadron z ~ "ri the squared hadron spatial transverse momentum with respect to the virtual 

photon direction p^ and the angle between the leptonic and hadronic planes 4> [l3| (see Fig. [5]) . Here Eq is the beam 
energy, E' and 9 are the scattered electron energy and angle, v = Eq — E' is the virtual photon energy in the lab 



4 



frame, and M is the proton mass. We will also use momentum transfer t = {q — Ph)'^, Feynman xp = 2p^'^' /W and 
the projection of the hadron momentum onto the photon direction p|| as alternative variables when they help with 

the physical interpretation. Here W = \/ {q + P)^ is the invariant mass of the final hadronic system and the CM 
label denotes the center-of-momentum frame. 




The CEBAF Large Acceptance Spectrometer (CLAS) in Hall B at Jefferson Lab allows us to study the five-fold 
differential semi-inclusive cross section over a large range of four-momentum transfer from 1.4 to 7 (GeV/c)^ and 
Bjorken x from 0.15 to 1 (see Fig. [31). CLAS enables us to measure distributions of the outgoing meson {z from 
0.07 to 1 and from 0.005 to 1.5 (GeV/c)^), in particular, full coverage in the azimuthal distributions that is very 
important. However, the covered kinematical interval is not rectangular in all five dimensions leading to a shrinkage 
of the four-dimensional acceptance when one of variables reaches its limits. The CLAS detector also has good particle 
identification capabilities, resulting in a clean selection of pious for this analysis. 



II. FORMALISM AND THEORETICAL EXPECTATIONS 



The unpolarized semi-inclusive cross section can be written in terms of four independent Lorentz- invariant structure 
fimctionsllal: 



dxdQ'^dzdp^c 



2t^o? Eh 



-K2(2 -y)J-cos -I- 2k cos 20^4 



(1) 



where inelasticity y — v/Eq, ^ 



2Mx 



, C = ^-y- jl'^y'^, e = = irpT? and Hi = Hi{x, z, Q'^^p^). In contrast to 



Ref. [Ill, we absorbed the \Jp\IQ'^ and p^lQ^ coefhcients in front of the structure functions H^, and Hi, respectively, 
into the structure function definition to let Hi reflect the recently identified leading twist contribution by Boer and 
Mulders p^ . Both structure functions include also an additional factor i to simplify relation with the azimuthal 
moments. 

In order to disentangle all four structure functions Hi, one has to measure the complete five- fold differential semi- 
inclusive cross section at a few different beam energies (as proposed in Ref. ^^20*1). In the present work we have data 
for only a single beam energy, and we relied upon the separation between the longitudinal (ctl ^ H2/ k ~ 2xHi) and 
transverse (ctt ~ 2xHi) cross sections performed in Ref. [2l[ and found to be compatible with R = 0.12 ± 0.06 (the 
weighted average over proton and deuteron data): 



Hi 1 (Tt 

H2 2kX CfT + 



1 



1 



2k.xI + R 



(2) 



where R = tl/ut is the longitudinal to transverse cross section ratio. 
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FIG. 3: Kinematical regions covered by the present experiment in different independent variables. Numbers on the plots give 
base-10 logarithm of the number of events. 



The azimuthal moments can be expressed in terms of the structure functions Hi as follows: 

(cose- - ''^ 



(3) 



(cos2(/)) = K- 



7^4 



712 + eHi 

These relations allow us to extract azimuthal moments from the data in order to determine Tis and H4. 
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The SIDIS cross section integrated over (j) and is given by [H, H^l : 



dxdQ^dz xQ^ 
Mxy 



1-y 



2E^ 



(4) 



where roman Hi are defined as caUigraphic Tii structure functions integrated over p^. In the parton model, the 
initial momentum of the struck quark is given by the proton momentum multiplied by the light-cone fraction x. If we 
instead consider the momentum carried by the struck quark after absorption of the virtual photon, then z represents 
the light-cone fraction of the momentum taken away by the produced hadron. In the region of forward going hadron 
(frame dependent) this cross section can be evaluated as the convolution of the parton density function f{x,Q^) 
obtained in inclusive processes and the parton fragmentation function D^{z, Q^) measured in e~^e~ collisions: 



H2{x, z, Q2) ^ e^^Mx, Q2) ® Dl^z, 



(5) 



where the sum runs over quark flavors i and is the charge of ith flavor quark. Instead, in the region of backward 
going hadron (frame dependent) the cross section is proportional to the fracture function M^{x,z,Q'^) uniquely 
defined in the semi-inclusive processes: 



H2{x, z, 



x{\-x)Ml'{x,z,Q'^ 



(6) 



The separation between the two processes is frame dependent and can only be studied phenomenologically. 

Values of the parton density function /(x, Q^) and the parton fragmentation function D^{z, Q^) can be obtained 
in pQCD inspired world data fits e.g. in Refs. [H, [2^ and [1^ [23, HI], respectively, but only Re f.[26l a llows for 
hadron charge separation. The fracture function is only studied for proton and neutron production [29l . |30| and tt"*" 
fracture function is completely unknown. 

In practice the measured cross section also depends on the transverse momentum pT of the hadron. The intrinsic 
motion of partons in the proton (Cahn effect [31]) leads to an exponential p^-behavior of the structure function 0,2'. 



7^2(0;, z,(5^py) 



i?2(x^,z) 
7r(p2,) 



■ exp 



~Pt/{pI 



The mean squared transverse momentum in the naive parton model is given by the sum of two terms [32, l33|, |34 



52 + a2z2 



(7) 



(8) 



where is the mean squared intrinsic transverse momentum of the partons, a^z^ is the mean squared parton 
transverse momentum transferred to the hadron, and is the mean squared transverse momentum acquired during 
fragmentation. 

The (/)-dependent terms in Eq. [1] also appear in NLO pQCD because the radiation of hard gluons leads to an 
azimuthal variation [2^ . However, this effect is expected to be important at energies higher than that of the present 
experiment [s^l • This is because the transverse momentum generated by the hard gluon is generally larger than that 
accessible in our experiment. In our energy range, the main contributions to the (/)-dependence of the cross sections 
are expected to be the Cahn and Berger effects for H3 and the Boer-Mulders function [lOj] for Ha (see also 
Refs. [H, [13, m, |3^ ) . The Cahn effect arises from the simple kinematics of partons with transverse momentum and 
can be calculated explicitly in the limits ^ 00 and z ^ 1 (see Refs. (s^ [ssl [s^. [40l|). 

The Berger effect is the exclusive production of a single pion from a free, struck quark that radiates a gluon, 
produces a qq pair, and recombines with the q. The formation of this pion through one-gluon exchange yields a cos 
dependence proportional to the hadron wave function. Since such production mechanism does not require any initial 
transverse momentum of struck parton it is completely orthogonal to the Cahn effect. 

Explicitly neglecting intrinsic parton transverse momentum one has |4H : 



Hz 



zhih - ^h) 



+ eHi 



rjll + iAz^ + ^^^y^lf 



(9) 
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and 

(11) 

with "0(0 being the pion wave function and rj = 1 + e/2x. 

The contribution of the Boer-Mulders function gives the probabihty to find a transversely polarized quark in the 
unpolarized proton. Explicitly in LO pQCD and p^/Q^ — > one has |19| : 

W4 _ l-y v\ ^Y..A^hi{x)Et\z) ^^^^ 



Here we defined 



and 



where ruh is the mass of the detected hadron, {x) is the momenuni distribution of transversely polarized quarks 
in the unpolarized proton (Boer-Mulders function) and H^^(z) is the Collins fragmentation function [i^l describing 
fragmentation of a transversely polarized quark into a polarized hadron. The Collins fragmentation function was 
parametrized using e+e~ data in Ref. \^^. 

These three main effects predict different kinematic dependencies. For example, the contribution of the Boer- 
Mulders function in 7^4 is of leading order, and therefore should scale with . On the other hand, both the Cahn 
and Berger effects have a non-perturbative origin and should decrease with rising Q^. The Cahn and Berger effects 
have opposite signs, but both increase in magnitude with z. The Berger effect should also increase in magnitude 
with X as the exclusive limit is approached, whereas the Cahn effect does not have any x dependence. To distinguish 
among these physical effects, one needs to perform a complete study of all kinematic dependencies in the data. 



III. DATA ANALYSIS 

The data were collected at Jefferson Lab in Hall B with CLAS [43| using a 0.354 g/cm^ liquid-hydrogen target and a 
5.75-GeV electron beam during the period October 2001 to January 2002. The average luminosity was 10'^^ cm~^s~^. 
CLAS is based on a six-sector torus magnet with its field pointing azimuthally around the beam direction. The torus 
polarity was set to bend negatively charged particles toward the beam line. The sectors delimited by the magnet coils 
are individually instrumented to form six independent magnetic spectrometers. The particle detection system includes 
drift chambers (DC) for track reconstruction [i^, scintillation counters (SC) for time-of-flight measurements [i^ . 
Cherenkov counters (CC) for electron identification and electromagnetic calorimeters (EC) for electron-pion 
separation [i^. The CLAS can detect and identify charged particles with momenta down to 0.2 GeV/c for polar 
angles between 8° and 142°, while the electron-pion separation is limited up to about 50° by the CC acceptance. 
The total angular acceptance for electrons is about 1.5 sr. The CLAS superconducting coils limit the acceptance for 
charged hadrons to about 80% at = 90° and about 50% at = 20° (forward angles). 

The electron momentum resolution is a function of the scattered electron angle and varies from 0.5% for < 30° up 
to 1-2% for 6 > 30°. The angular resolution is approximately constant, approaching 1 mrad for polar and 4 mrad for 
azimuthal angles. Therefore, the momentum transfer resolution ranges from 0.2 to 0.5 %. For the present experiment 
the invariant mass of the struck proton (W = \/ {P + q)^) has an estimated resolution of 2.5 MeV for beam energies 
less than 3 GeV and about 7 MeV for larger energies. In order to study all possible multi-particle states, we set the 
data acquisition trigger to require at least one electron candidate in any of the sectors, where an electron candidate 
was defined as the coincidence of a signal in the EC and Cherenkov modules for any one of the sectors. 



A. Generic procedures 



Both the e and 7r+ were detected within the volume defined by fiducial cuts. These geometrical cuts selected regions 
of uniform high efficiency by removing areas near the detector boundaries and regions corresponding to problematic 
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SC counters or DC readout. For electrons the fiducial volume limitations are mostly due to the Cherenkov counter, 
which is necessary for electron identification, and the electromagnetic calorimeter, which is used in the trigger. The 
CLAS Cherenkov counter's optics reduce significantly its azimuthal acceptance, in particular in the region of small 
polar scattering angles, where the light collection mirrors are small. Moreover, the Cherenkov counter extends only up 
to ~50° in the polar scattering angle of an inbending charged particle. The trigger threshold for the electromagnetic 
calorimeter limits the lowest electron momentum, which in our case was about 0.64 GeV/c. 

CLAS achieves its best charged-particle acceptance for 7r+, since complete identification requires only information 
from the drift chambers and the scintillation counters, which are limited in coverage only by the CLAS torus magnet's 
coils. For the standard torus configuration, 7r+ particles bend outward toward larger angles, where the useful detector 
area between the coils is greater. 

Small corrections to the momenta of the and 7r+ were necessary because of distortions in the drift chambers 
and magnetic fields not accounted for in the tracking routines. Correction parameters were determined by minimizing 
the difference in the missing mass for ep ep and ep eTT~^n from known values (see Ref. [49|). The magnitude of 
these kinematic corrections was well below the CLAS resolution leading to sub-percent changes in the measured cross 
section. 

Events were selected by a coincidence of an electron and a Tr"*" whose identification criteria are described in the next 
section. The trigger gate time in CLAS was 150 ns, but, due to the limited range of particle momenta, the effective 
time window for a coincidence was much smaller. This, and the relatively low beam current in CLAS (about 7 nA), 
ensured that accidental coincidences were negligible. 



The electron identification is based on combined information from the CC, EC, DC and SC. The fastest (as measured 
by the SC) negatively charged (as determined from DC tracking) particle having EC and CC hits is assumed to be 
an electron. However, the large rate of negatively charged pions contaminates the sample of reconstructed electrons, 
in particular in the region of low momenta and large polar scattering angles. Moreover at lowest accessible polar 
scattering angles CC efficiency is reduced due to geometrical constraints. This contamination can be eliminate d by 
using SC and DC information to better correlate the particle track and the time of the SC hit with the CC signal 5(jl|. 




We estimated the electron efficiency after this process to be greater than 97% and the corresponding inefhciency was 
propagated to e~ identification systematic uncertainty. The inefficiency is maximum at the lowest Q^. 

The CC becomes less efficient at distinguishing electrons from pions for momenta above the Cherenkov light 
threshold for pions (|p| ~ 2.7 GeV/c). However, in this kinematic region the EC signal can be used to remove the 
remaining pion contamination. The minimum- ionizing pion releases a nearly constant energy in the EC, independent 
of its momentum, whereas an electron releases an almost constant energy fraction of about 30% in the EC. Fig. 3] 
shows a contour plot of events with momentum p determined from the DC and total energy in the EC normalized by 
|p|. Pion-electron separation, in this case, increases with particle momentum. 

Pion identification is based on time of flight as measured with the SC and momentum as measured with the 
DC. Since the distance between the target and SC is independent of the scattering angle, the efficiency of pion 
identification depends only on the pion momentum and therefore on z. The time-of- flight resolution decreases with 
pion momentum leading to larger pion identification inefficiency. A contribution proportional to this inefficiency was 
added to 7r+ identification systematic uncertainty. Furthermore, the time-of-flight interval between different hadron 
species decreases with hadron momentum resulting in larger contamination. 

Fig. [5] shows how effectively this procedure removes the background under the exclusive n^n peak without a 
significant loss of good events. The example of the exclusive n~^n peak is important because these pions have large 
momenta, which makes their separation by time of flight more difficult than in the semi-inclusive case, where the 
pions have slightly lower momenta. 

A positively charged particle identified as a pion may in some cases be a positron from e^e~ pair production. This 
background becomes important at low momenta and at </) « or </) « 180°. To remove this contamination we applied 
the cut 



where Mtof is the mass of the positive particle as measured by the TOF and M(e~/i+) is the invariant mass of the 
measured system of two particles in GeV/c^ (assuming them to be e+ and e~ and gtof = 0.01 (GeV/c^)^). The 
remaining contribution from e+e~ pairs is negligible over the entire kinematic range after the cut. 

The electron, detected in coincidence with the pion, may be a secondary electron, whereas the scattered electron is 
not observed. To remove this background we measured lepton charge symmetric e"'"7r+ coincidence cross section from 



B. Particle Identification and Backgrounds 




M^{e-h+) > 0.012 exp -M^of/C^'^tof) 



(14) 
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FIG. 4: Contour plot of particle momentum p from tracking versus particle energy deposited in the calorimeter ECtot normalized 
by IpI . Events on the left correspond to pions and those on the right to electrons. Only fiducial cuts were applied. The Cherenkov 
detector, providing the basic electron identification, allows to identify clearly the electrons up to |p| « 2.7 GeV/c. The dashed 
line shows the cut applied to the data to remove remaining pion contamination at |p| > 2.7 GeV/c. 




(GeV/c ) 



0.6 0.8 1 1.2 1.4 
(GeV/c^)^ 



FIG. 5: Measured squared mass of positive hadrons (left) and squared missing mass for the e /i"'"-system in the region of the 
tt"*"?! exclusive peak (right). The shaded area indicates hadrons identified as pions. The exclusive Tr+n was removed from the 
analysis by a cut Mi > 1.08 {GeV/c^f. 



the same data and removed its contribution from our data. This contamination is hmited to a few lowest-x points 

where it achieves 5% at most. 

Another source of contamination comes from production at high hadron momenta. At low hadron momenta the 
TOF system is able to distinguish pions from kaons, but above |p;j| ~ 1.2 GeV/c the peaks of the two particles begin 
to mix. However, large hadron momenta make two-kaon production less likely due to the correspondingly high energy 
threshold, and therefore most of the background comes from single kaons associated with A and TP production. In 
order to suppress the kaon contamination we applied two cuts: a kinematical cut that removes A, TP and A*(1520), 
and a TOF cut < m'^ + 2aM^(TOF) that suppresses low- momentum kaons. The mass resolution of the TOF system 
was determined by fitting the width of the pion peak, which yielded crM^(TOF) = 0.022|p^| exp (0.6-^/|pJ), where is 
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given in GeV/c and <tm^(tof) in (GeV/c^)^. Corrections for the remaining kaons from semi-inclusive production above 
the two kaon threshold were made using the ratio of to 7r+ semi-inclusive cross sections obtained from a pQCD- 
based Monte Carlo (MC) event generator (see the following section), weighted with the kaon/pion rejection factor 
obtained from the simulation itself. Kaons from the MC were propagated through the entire chain of the reconstruction 
procedure exactly in the same way as was done for pions, and the fraction f{K^) of kaons reconstructed as pions 
was obtained. This number was normalized to the fraction /(tt^) of simulated pions reconstructed by the procedure. 
This kaon/pion rejection factor was parameterized as a function of the hadron momentum. The contribution from 
the background varied from to 20% with an average of 1%, and our procedure reduced the kaons by a factor of 
two at 2.3 GeV/c with an increasing reduction factor at lower hadron momenta. 

C. Empty target contribution 

Empty target runs were analyzed in exactly in the same way as the full target runs and subtracted from the data to 
eliminate scattering from the target endcaps. The total charge collected on the empty target is an order of magnitude 
smaller than the one for the full target data. In order to increase the statistics of the empty target distributions, 
we made the assumption that the ratio of full to empty target event rates factorizes as a function of all variables. 
Thus one can obtain the ratio of empty to full target rates (ranging from to 18% with an average value of 4.7%) for 
the five-fold differential cross section as a product of one-fold differential ratios. The contribution of empty target is 
typically smaller than the total systematic uncertainty, reaching its maximum at the two pion threshold. 

D. Monte Carlo Simulations 

Detector efficiencies and acceptances were studied with a standard CLAS simulation package GSIM [51]. The 
simulated data obtained from GSIM can then be analyzed using the event reconstruction routine exactly in the same 
way as the measured data. This allows a complete determination of the detector efficiency plus acceptance. 

The first step of the simulation is to generate e~7r+ coincidence events based on a pQCD-like SIDIS parameteri- 
zation (5^ at leading order for the semi-inclusive contribution and on the MAID2003 model [53| extrapolated to the 
W > 2 GeV/c^ region with the parameterization from Ref. [S^l for the exclusive 7r+n reaction. Distributions of counts 
from the experimental data and GSIM simulations are shown in Fig. [6l The same cuts are applied to both data and 
MC as described in the previous section. 

The MC yield reproduces the shape of the experimental data fairly well. In order to keep systematic uncertainties on 
the acceptance plus efficiency small (we estimated them to be 10%) we had to extract fully-differential cross sections 
in narrow kinematic bins. Bins with combined acceptance and efficiency < 0.1 %, corresponding to 25% filling per 
each dimension, were discarded. The average value of combined acceptance and efficiency was about 25%. 

To test Monte Carlo simulations of electrons we extracted the inclusive structure function F2 and compared it to 
the world data in our kinematic range. An example of this comparison at = 2 (GeV/c)^ is shown in Fig. [T] 

For positively charged hadrons we tested Monte Carlo simulations by extracting the elastic scattering cross section 
measured in electron-proton coincidences. The normalized event yield was compared to normalized GSIM simulation 
yield, based on form-factors from Ref. (56j . The obtained ratio, shown in Fig. [HI is in good agreement with unity in 
the central region of Q^, but rises at large due to unresolved inelastic contamination. 

Furthermore, the efhciency of 7r+ production reconstruction in the present data set was tested in the measurement 
of the exclusive pion production published in Ref. 49]. 

E. Binning 

The data were divided into kinematic bins as follows: 

• g2 - 10 logarithmic bins (with centers at) 1.31-1.56 (1.49), 1.56-1.87 (1.74), 1.87-2.23 (2.02), 2.23-2.66 (2.37), 
2.66-3.17 (2.93), 3.17-3.79 (3.42), 3.79-4.52 (4.1), 4.52-5.4 (4.85), 5.4-6.45 (5.72), 6.45-7.7 (6.61) (GeV/c)^; 

• X - 25 logarithmic bins in the interval from 0.01 to 1; 

• z - 25 logarithmic bins in the interval from 0.01 to 1; 



• PT - 10 logarithmic bins (with centers at): 0-0.1 (0.07), 0.1-0.2 (0.16), 0.2-0.3 (0.26), 0.3-0.41 (0.36), 0.41-0.53 
(0.47), 0.53-0.65 (0.58), 0.65-0.8 (0.71), 0.8-0.97 (0.86), 0.97-1.22 (1.04), 1.22-1.64 (1.25) GeV/c; 
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FIG. 6: Comparison of e tt"*" coincidence data (full triangles) and MC raw yields (open circles) as a function of one of the 
kinematic variables. The other variables were kept fixed at = 2.4 (GeV/c)^, x = 0.30, z = 0.37, pr = 0.22 {GeV/cf. The 
MC simulation yield was normalized to the integrated luminosity of the experiment. Error bars are statistical only and they 
are smaller than the symbol size. 



• (j) - 18 linear bins in the interval from to 360°. 

The bin sizes were chosen large enough to accommodate CLAS angular and momentum resolutions reducing bin 
migrations, but they were small enough to avoid reconstruction efficiency model dependence. Indeed the average bin 
migration is about 30%, consistent with expectation for bin size close to la of detector resolution. The centers of the 
and pr bins coincide with the mean values of these variables in the raw data after all cuts but before acceptance 
corrections. 

The general rectangular grid described above was used to sort the data. But the measured kinematic volume is not 
rectangular. Hence a fraction of bins in one dimension can result empty depending on the values of the other four 
variables. This is in particular the case when one of variables is close to the kinematic limits of the accessible region. 
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FIG. 7: Inclusive structure function F2 ai — 2 (GeV/c)'^ extracted from the present experiment (full triangles) in comparison 
with previous world data (open circles) from Ref. [S^ and references therein. The curve is a parameterization from Ref. [55|. 
The error bars are statistical only and an overall systematic uncertainty for the present experiment of the order or 5% is 
estimated. 
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FIG. 8: Ratio of e p coincidence events from the data and GSIM Monte Carlo simulations. The simulation event generator 
was based on proton form- factors from Ref. [5^. The coincidences were selected by the following set of cuts: \W'^ — M^\ < 0.2 
{GeV/c^f, \Mj,\ < 0.01 (GeV/cy and \\(fih - (l>e\ - n\ < 4 de grees. The error bars are statistical only and an overall systematic 
uncertainty of the order or 10% is estimated. Error bars are smaller than the symbol size. 



F. Five- fold differential Cross Section 



The five-fold differential semi-inclusive cross section was extracted for each kinematic bin from the number of 
measured events Ndat according to the relation: 



d^a G. 



dat 



dxdQ'^dzdp'^dct) AxAQ^AzAp'^A(j> 

Feff/acc{x,Q'^,Z,p'^,(l)) 



(15) 
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with data inverse luminosity given by 



Gdat - 7 - N TD ' ^^^^ 



p — 0.0708 g/cm^ is the hquid- hydrogen target density, L = 5 cm is the target length and Qfc is the total charge 
collected in the Faraday Cup (FC), corrected for dead time. Ff.fffacc{x,Q'^,z,p^,(j)) is the acceptance/efficiency 
correction obtained with Monte Carlo simulations: 



Nrecix,Q^,Z,p'^,(f)) 
(JM{x,Q'^,Z,p\,,(f>) 

with its normalization factor given by 

Gsim = ■ (18) 

Here Nrec{x,Q'^ , z,p'^,(l)) is the number of Monte Carlo events reconstructed in the current bin, (74/(2;, Q^, z,^!^, 0) 
is the cross section model used in the event generator, r is the complete phase space volume of event generation and 
Ntot is the total number of generated events. 

The final cross sections were corrected for radiative effects using the analytic calculations described in Refs. [s^-ISTj 
implemented in the Monte Carlo generator. It includes both radiative corrections to the SIDIS spectrum and the 
radiative tail from exclusive TT~^n production. The average contribution of radiative corrections is about 6% with 
largest contribution close to the two pion threshold. The magnitude of radiative corrections increases with z and p^. 



G. Azimuthal dependence 

A separation of the constant, cos cj) and cos 20 terms in Eq. [T]has been performed using two methods, either a fit to 
the ^-distributions or an event-by-event determination of azimuthal moments. Both methods should give compatible 
results. By studying the two methods in detail we concluded that both give unreliable results if the (/)-distribution 
contains regions of poor detector acceptance. Therefore we excluded kinematic points where the (/)-acceptance was 
inadequate. This reduced significantly the kinematic range of the extracted moments. Nevertheless, the kinematic 
bins with incomplete (/)-coverage can still be used in a multidimensional fit exploiting continuity in the other variables. 

In the first method we fit the ^-distribution to the function (To(1 + 23 cos (f) + 2Ccos20) using MINUIT [s^ and 
extracted the coefficients ctq, B and C and their statistical uncertainties. These coefficients give the ^-integrated cross 
section, (cost/)) and (cos2(/)), respectively. 

The second method of moments was used in a previous CLAS paper [5^ , but due to the strong effect of acceptance 
on even moments, we developed the necessary corrections described below. The fully differential cross section can be 
written as: 

CT = + cos + V2 cos 2(t> . (19) 
The acceptance/efficiency correction can be expanded in a Fourier series in (f> as: 



oc 



Fef f/acc = + X! ^" '^'^^ n(j) + '^Bn s'm ncj) . (20) 



2 

n—l n—1 



The coefficients i?„ are fairly small in CLAS. Fourier spectrum of the raw data and MC yields is shown in Fig. O 
Let us define the Fourier coefficients of raw yield (before acceptance/efficiency correction) 



1 /•^'^ 

y„ = — / a {4>) cos n(j)d(j) , (21) 

where a = F^ff/acc'^ is the cross section distorted by the acceptance/efficiency correction F^f f/acc- Then combining 
Egs. fTOll^ andl^we obtain an infinite series of hnear equations relating Fourier coefficients of the raw yield and the 
physical cross section: 

r„ = A^V, + ^Ilzl±Al±lv, + ^^l^l±All±lv, , (22) 
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FIG. 9: Extracted Fourier components of the raw data (full triangles) and the Monte Carlo (open circles) yields. Error bars 
are statistical only and they are smaller than the symbol size. 



where An — v4_„ and n = 0, 1, 2, . . . . The magnitudes of An and Yn decrease rapidly with n and are consistent with 
zero for n > 10 (see Fig. As one can see in the Figure the Monte Carlo simulations describe fairly well Fourier 
spectrum of the data, in particular for large n. Therefore, we can cut the infinite set of equations for MC Yn at some 
arbitrary n = N and solve the resulting system of N linear equations to obtain An coefficients for n = 0, 1, 2, ... TV. 
Assuming that GSIM reproduces the CLAS acceptance and efficiency (within systematic uncertainties treated later), 
coefficients An should be the same as in the data. We used these efficiency /acceptance Fourier coefficients A„ in 
the expression for data Yn to extract the measured cross section (/)-terms: Vq, Vi and V2. We fit the overdetermined 
system of N linear equations with these 3 unknowns using the weighted linear least squares fitting routine TLS in the 
CERNLIB library 

The stability of the solution as a function of N is shown in Fig. [TOl From this plot we concluded that = 7 is the 
minimum number of moments necessary to extract sensible (cos 0) and (cos 2(p) for the present kinematics. In the 
following we made the more conservative choice of A^ = 20. 

A typical acceptance-corrected (/)-distribution is shown in Fig. [11] together with the two methods of extracting 
moments, which are in good agreement. The systematic uncertainties on the ^-dependent terms are larger than the 
difference between the two methods (see the following section). 

The acceptance correction was also tested by dividing each bin in two parts by cutting the corresponding scattered 
electron energy range in two equal intervals and comparing the extracted (cos and (cos 2(/)) terms in these two 
acceptance regions. An example of this test is shown in Figs. [T2] and fT3l It can be see that the acceptance corrections 
are significant and sometimes differ for the two separated regions, however the final reconstructed (cos 0) and (cos 20) 
comes out to be consistent within statistical and systematic uncertainties. 

H. Systematic Uncertainties 

The systematic uncertainties for the measured absolute cross sections are considerably different from those for 
the azimuthal moments, because many quantities drop out in the ratios measured by moments. Most of systematic 
uncertainties are point-to-point correlated and evaluated on a bin-by-bin basis with the exception of the overall 
normalization, efficiency, and radiative and bin-centering corrections, for which a uniform relative uncertainty was 
assumed. The following sections discuss these uncertainties. 
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FIG. 10: Stability of the V2 cross section term as a function of the number of moments taken into account in the extraction 
procedure. Error bars are statistical only. 




100 150 200 250 300 350 
degrees 



FIG. 11: The (^-dependence of the data taken at = 2 {GeY/cf, x = 0.24, 2 = 0.18 and Pr = 0.5 {GeY/cf (full triangles) 
together with the results of the azimuthal moment (solid lines) and fitting (dashed line) methods. Error bars are statistical 
only. 



1. Cross Section 



The total systematic uncertainties on the five-fold differential cross section vary from 11 to 44% with a mean value 
of 16%. Apart from systematics due to the efficiency corrections discussed in Section fill Dl the major contributions 
come from detector acceptance and electron identification. The relative value of most uncertainties is amplified at 
the two pion threshold, where the cross section vanishes. 

The estimate of the systematic uncertainty from efficiency modeling comes from a comparison of cross section 
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FIG. 12: The p|^-dependence of the {cos<j!>) taken at = 2 (GeV/c)^, a; = 0.30 and z = 0.21 obtained from two different 
detector acceptance regions (triangles and squares). Full markers correspond to the data before the acceptance correction 
and the empty markers show acceptance corrected (cos(/>). The two data sets are shifted equally along the x-axis in opposite 
directions from their central values for visibility. Error bars are statistical only. 
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FIG. 13: Same as Fig. [T2l except for (cos2(^). 



extractions using two different event generators: one uses a LO pQCD model, while the other is based on the sum 
over several exclusive channels. 

The systematic uncertainties on the acceptance were estimated from the variation in the absolute cross sections 
obtained using each of six CLAS sectors separately to detect the electron (pion) and then integrating over the pion 
(electron) wherever else it appeared. This uncertainty was estimated bin-by-bin and reflects the ability of Monte 
Carlo to describe the detector non-uniformities. The uncertainty increases at low polar scattering angle, and therefore 
low-Q^ for electrons and low-p^ for pions, where the azimuthal acceptance of CLAS is reduced. 
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Systematic uncertainties arising from electron identification were estimated by comparing two different methods 
(as in Ref. [13] ) of pion rejection, one based on Poisson shapes of Cherenkov counter spectra and another on the 
geometrical and temporal matching between the measured track and Cherenkov signal. This uncertainty appear 
mostly at low-Q^, where CC is less efficient. 

The systematic uncertainty arising from 7r+ identification has two contributions. One was estimated from the 
difference between the ratios of events in the missing neutron peak before and after pion identification as calculated 
for data and GSIM simulations. The second part comes from our treatment of kaon contamination which was assumed 
to be 20%. The two errors were added in quadrature. 

Radiative corrections are model-dependent. To estimate this systematic uncertainty we changed the model used in 
the radiative correction code by 15% and took the resulting difference as an estimate of the uncertainty. 

There is an additional overall systematic uncertainty of 1% due to uncertainties in the target length and density. 
The target length was 5±0.05 cm and the liquid-hydrogen density was p = 0.0708±0.0003 g/cw? giving approximately 
a 1% uncertainty. 

The systematic uncertainty on the bin centering correction was estimated in the same way as for the radiative 
corrections. The model was changed as described above and the difference between the two centering corrections was 
taken as the uncertainty. 

The empty target subtraction introduces a small systematic uncertainty due to the assumption of cross-section- 
ratio (empty to full target) factorization in the individual kinematic variables. This uncertainty was estimated by 
comparing the factorized and direct bin- by-bin subtraction methods. 

These main contributions are listed in Tabic |T1 All systematic uncertainties shown in the table were combined in 
quadrature. 

TABLE I: Systematic uncertainties of the semi-inclusive cross section. 



Source 


Variation range 

% 


Mean value 

% 


Overall normalization 


1 


1 


e~ identification 


1.8-13 


3.5 


tt"*" identification 


0.9-6.7 


2.1 


e~ acceptance 


0-19 


5.3 


TT^ acceptance 


0-52 


4 


Efiiciency 


10 


10 


Radiative corrections 


2 


2 


Empty target subtraction 


0-0.7 


0.2 


Bin centering correction 


0.7 


0.7 


Total 


11-54 


14 



2. Azimuthal Moments 

Azimuthal moments (see Eq. [3|) have the advantage of smaller systematic uncertainties since many of them cancel in 
the ratio. In particular, systematic uncertainties of overall normalization, kinematic corrections, particle identification, 
efficiency, empty target subtraction and bin centering cancel. The remaining systematic uncertainties are due to non- 
uniformities in the CLAS acceptance and radiative corrections. The uncertainties due to CLAS acceptance were 
estimated as the spread between the central values of the azimuthal moments obtained using each single CLAS sector 
to detect the electron or pion and then integrating over the second particle (pion and electron respectively). This way 
we obtained the influence of the electron and pion acceptances se par ately. Similar conclusions about the acceptance 
influence on the azimuthal moment extraction were made in Ref. [63| . 

To estimate the systematic uncertainties of the radiative corrections, we made a few calculations in randomly 
chosen kinematic points comparing correction factors obtained with our model, changing by 15% the exclusive 7r+ri 
contribution or modifying by 30% the 7^3 and 7^4 structure functions. The difference in the correction factor was 
taken as the estimate of this systematic uncertainty. The variation range and averaged value of these systematic 
uncertainties are given in Tables [ill and IIIII for the (cost/)) and (cos2(/)) moments, respectively. 
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TABLE II: Systematic uncertainties of {cos(j)). 



Source 


Variation range 


Mean value 


e~ acceptance 


0-0.06 


0.016 


7r+ acceptance 


0-0.13 


0.016 


Radiative corrections 


0.005 


0.005 


Total 


0.005-0.13 


0.026 



TABLE III: Systematic uncertainties of (cos 20). 



Source 


Variation range 


Mean value 


e~ acceptance 


0-0.08 


0.015 


7r+ acceptance 


0-0.12 


0.011 


Radiative corrections 


0.003 


0.003 


Total 


0.003-0.12 


0.021 



3. Structure Functions 

One additional systematic uncertainty appears in the extraction of the structure function 7^2 from the measured 
combination H2 + eTi-i. In this case some transverse to longitudinal cross section ratio R should be assumed. In our 
results on the structure function 7^2 we included a 50% systematic uncertainty on R. This does not affect strongly the 
extracted structure function 7^2 (see Eq. [2]), in the same way as the inclusive structure function F2 is weakly sensitive 
to the ratio R for forward-angle scattering. The assumed 50% precision leads to the systematic uncertainty shown in 
the Table HVl 

TABLE IV: Additional systematic uncertainty on Ti.2- 



Source 


Variation range 

% 


Mean value 

% 


R ratio 


0.6-1.9 


1.5 



IV. RESULTS 

The obtained data allow us to perform studies in four different areas: hadron transverse momentum distributions, 
comparison of the ^-independent term with pQCD calculations, search for the target fragmentation contribution and 
study of azimuthal moments. We present these analyses in the following sections. 

A. Transverse Momentum Distributions 

The ^-independent part of the cross section falls off exponentially in as shown in Fig. [TH This has been 
predicted in Ref. [sij to arise from the intrinsic transverse momentum of partons. We observe no deviation from this 
exponential behavior over the entire kinematic domain of our data. 

By studying the p|,-dependence in our data at various values of z, we have extracted the z-dependence of the mean 
transverse momentum (p^), defined within the Gaussian model, in Eq. [71 and obtained by fitting ^^-distributions in 
each (cc, Q^, z) bin. Fig. [TSl shows a clear rise of (p^) with z. We compared this with the distribution given in Eq. [8] 
with a? = 0.25 and &^ = 0.20 (GeV/c)^ based on previous data [H, [s^lsll. Significant deviations from this behavior 
were found at low-z, which can be explained as a threshold kinematic effect. The maximum achievable transverse 
momentum p™"^ ~ zv becomes smaller at low z, because v is limited by the 5.75-GeV beam energy, and p™"^ is 
smaller than the intrinsic transverse momentum of partons which is at first order independent of beam energy. This 
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FIG. 14: The py-dependence of the (^-independent term TI2 + tHi at s = 0.24 and z — 0.30. The hnes represent exponential fits 
to the data for = 1.74 (GeV/c)^ (full circles and solid line), = 2 (GeV/c)^ (full squares and dashed line), and — 2.37 
(GeV/c)'^ (triangles and dotted line). The errors bars are statistical only. 



leads to a cut on the py-distribution, which is not present in high energy experiments. To account for this low-energy 
effect we modified the parameterization as: 

The dotted curve in Fig. [13] shows that this new parameterization follows the data points, but the absolute normal- 
ization given by the parameters a and 6 is still too high. This modification breaks the factorization between x, 
and pT in the low-z region because the p|n-distribution now depends also on x and . 




FIG. 15: The z-dependence of (py) at = 2.37 (GeV/c)^ and x — 0.27. The points are the data from the present analysis. 
The curves show the maximum allowed = (p^^)"""^ (dashed), the parameterization of high energy data from Eq. [8] (solid), 
and the low-2 modification from Eq. [531 The error bars are statistical only and they are smaller than the symbol size. 

At large z, is also large. Therefore, we can check the factorization of from x and Q^. Fig. [TBI shows 
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no appreciable dependence of the mean transverse momentum for x < 0.5 corresponding to the missing mass 
M|. < 1.6 (GeV/c^)^, i.e. the A resonance region. 
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FIG. 16: The x-dependence of (py) at — 2.37 (GeV/c)^ and z — 0.34. The points are from the present analysis. The curves 
show = (p|h)'""^ (dashed) and a constant fit to the data (solid). The error bars are statistical only. 

The transverse momentum distribution exhibits a small variation with over the covered kinematic interval as 
seen in the different slopes in Fig. [TH However, the coverage is insufficient to observe the logarithmic pQCD 
evolution of (p^) with discussed in Ref. [6l|. 

B. Comparison with pQCD 

In order to compare the (/)-independent term with pQCD predictions, we assumed a constant longitudinal to trans- 
verse cross section ratio i? = 0.12 (2l| . 

Since there is no TMD-based approach to which we could directly compare our data, we integrated the measured 
structure functions 7^2 in in order to compare Ti.2 measured in this experiment with H2 from pQCD calculations. 
We integrated Eq. [1] in and and compared with Eq. [5] obtaining 



(Pt 



,Pt) 







(24) 



defined by the pion threshold, which limits the longitudinal hadron momentum in the lab frame to 



where the upper limit of integration is given by the smaller of the quantities (p|,)™°^ = (zv)'^ — m\ and the value 



>pr-^^\[^^^n-M^)+Q^- 

2Mv{\ ~ z)-ml + 2zv'^ - 2Af„TO^| 



(25) 



(pp™"'. If we exploit the exponential behavior of the measured structure function 7^2 in 



This limits pf, < |p^p 
(see Eq. [T]) , the integration can be performed analytically leading to 



H2{x,Q'',z)^V{x,Q'',z)Ehe ^ 



(26) 
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where V{x,Q'^ , z) is the pT-independent part of the structure function and Erfi is the imaginary error function. By 
neglecting the factor in Eq. [l]and by extending the integral to infinity (as typically done in SIDIS analyses, 

see Eq. ^ , we find 

H2{x,Q\z)^V{x,Q^,z) . (27) 
In Figs. [TTI [TSl fT9l and [20l our integrated structure function H2 is compared to pQCD calculations given by: 

^^2(x,Q^z)- (28) 

where crj^"*^^^ is the hard scattering cross section for incoming parton i and outgoing parton j given in Ref . j62i] , is 
the parton distribution function for parton i taken from Ref. [23j , D'^^ is the fragmentation function for parton j and 
hadron 7r+ taken from Ref. \20\, and /i is the factorization/renormalization scale. These next-to-leading order (NLO) 
calculations include a systematic uncertainty due to arbitrary factorization/renormalization scale variations [63}, 
indicating the size of possible higher order effects. This was evaluated by variation of each scale by a factor of two 
in both directions and the obtained differences for all scales were summed in quadrature. NLO calculations within 
their uncertainty lie closer to the data in the low-z region than leading order (LO) ones. The difference between the 
data and NLO pQCD is at most about 20%. At low x and z < 0.4 the data are higher than NLO calculations, while 
at largest x both the LO and NLO calculations lie above the data. The multiplicity ratio H2/F2 shown in Fig. [21] 
demonstrates the same level of agreement between data and pQCD calculations as H2 alone. This suggests that the 
differences between the data and theory do not cancel in the ratio. 

The widening systematic uncertainty band in the NLO calculations at high x suggests that the discrepancy with 
the data here might be due to a significant contribution from multiple soft gluon emission, which can be resummed to 
all orders in as in Refs. [g^, [fi^- Similar results were obtained in Ref. ^66] from the comparison between HERMES 
7r+ SIDIS data and NLO calculations. 

The difference between the data and calculations in some kinematic regions leaves room for an additional contribu- 
tion from target fragmentation of <20%. However, the possible presence of higher twists at our relatively small 
values casts doubt on the attribution of data/pQCD differences to target fragmentation. In order to better explore 
target fragmentation, we studied the t and Xi^-dependencies of Ti.2 as described in the following section. 

The pQCD calculations are significantly biased by the assumption of favored fragmentation [2^. In fact, using 
unseparated -\- h~ fragmentation functions as directly measured in e^e~ collisions, one obtains curves that are 
systematically higher by about 20%. 

C. Target Fragmentation 

In Leading Order (LO) pQCD, the structure function H2 is given by: 

H2{x, z, Q2) ^ elx[f,{x, Q^)D'^{z, Q^)A{0^h = 0) (29) 

i 

+ (1 - x)Ml\x, z, Q^)Ai9^h - 7t)] , 

where D'^{z) is the fragmentation function, M^{x, z) is the fracture function 1] and A{9^h) is the angular distribution 
of the observed hadron [bJ- The fracture function M^{x,z) gives the combined probability of striking a parton of 
flavor i at X and producing a hadron h at z from the proton remnant. This function obeys the pQCD evolution 
equations [H, [131 similar to those for fi{x) and D^{z). The factorization of the hard photon-parton scattering and a 
soft part described by M^{x, z) has been proved in Refs. [1, 

Because the agreement between pQCD calculations and our data, shown in Fig. [171 was rather poor we could 
explore only qualitative behavior of the structure functions to search for the target fragmentation contribution. 

To estimate target fragmentation we used two alternative sets of variables: 1) z and t, where the squared 4- 
momentum transfer t provides added information on the direction of and 2) xp and which included the sign 
of the longitudinal hadron momentum in the center-of- momentum (CM) frame through Feynman xp- 

Target fragmentation is expected to appear at small z, where hadrons are kinematically allowed in the direction 
opposite to that of the virtual photon. In analogy with vector meson photoproduction measurements, a contribution 
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FIG. 17: The ^-dependence of H2 at = 2.37 (GeV/c)^ The data are shown by fuU triangles. The error bars give statistical 
and systematic uncertainties combined in quadrature. The solid line shows LO pQCD calculations using the prescription from 
Ref. [gJ, CTEQ 5 parton distribution functions 23], and the Kretzer fragmentation functions [26|. NLO calculations performed 
within the same framework (using CTEQ 5M PDFs) are shown by the shaded area, for which the width indicates systematic 
uncertainties due to factorization and renormalization scale variations 16311. 




FIG. 19: Same as Fig. 1171 except with H2 plotted as a function of rather than z at a; = 0.33. 
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FIG. 21: Same as Fig. 1171 except with H2/F2, where F2 is the inclusive structure function obtained in the same experiment. 
And NLO calculations are shown by the dashed line, without systematic uncertainties due to factorization and renormalization 
scales. 
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FIG. 22: Schematic representation of the leading particle target fragmentation mechanism. The hatched blob represents the 
structure function of a particle or a set of particles exchanged. In case of production the first particle in the blob is the 
neutron. 

from target fragmentation may come from w-channel exchange [6§| by a particle or a set of particles (see Fig. [22]) . In 
this case the cross section would be proportional to the structure function of the exchanged particle (e.g. a neutron) 
or set of particles [1^. In this case one would expect a peak at |t| = |i|maa;, m addition to the dominant peak at 
\t\ = \t\min due to Rcggc exchange in the t-channel. This u-channel production can be called the "leading particle" 
contribution in the target fragmentation region because the produced hadron carries almost all of the spectator 
momentum. However, the measured i-distribution shown in Fig. [23] displays the exponential behavior expected in 
Regge theory but does not show any evidence of the second peak at \t\ = \t\max- In Fig. [231 the solid line shows an 
expected it-channel exchange contribution assumed to be 1% of the t-channel term. As one can see this assumption 
is not supported by the data. This observation is in agreement with a known phenomenological rule that a particle 
not present in the initial state cannot be the leading particle in a target jet |7Q] . 
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It I (GeV/c)^ 

FIG. 23: The t-dependence of the (^-independent term 7^2 -I- eWi at = 2 {GeV/cf, x = 0.24 and z = 0.18. The data are 
shown as solid squares. The curves represent the exponential fit to the data (dashes) and the expected behavior of leading- 
particle target fragmentation (solid), assuming it to be 1% of the t-channel exchange term. The error bars are statistical only 
and they are smaller than symbol size. 

Another contribution may come from soft fragmentation of the spectator diquark. One can naively define all 
hadrons produced in the direction of the struck quark to be in the current fragmentation region, whereas those 
produced in the direction of the spectator diquark to be in the target fragmentation region. Since this definition is 
clearly frame-dependent, in the following we will use the CM frame. 

Fig. \7M shows the data for four pT bins as a function of a;^. They exhibit a wide distribution centered at xp — 0, 
which corresponds to the center of momentum. Such behavior is in good agreement with that observed in semi- 
inclusive 7r+ production by a muon beam at much higher energies [7l|. According to our definition, all hadrons at 
xp > come from current fragmentation, while those at xp < come from target fragmentation. 

In the CM frame, z mixes backward-angle production with the production of low-momentum forward-going 
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FIG. 24: The 2;F-dependence of the ^-independent term H2 + eHi at = 2 (GeV/c)^ and x — 0.26. The data are compared 
to LO pQCD calculations combined with a Gaussian pT-dependence from Eq. [7] for = 0.005 (GeV/c)^ (full circles and 
solid curve), = 0.13 (GeV/c)^ (full triangles and dashed curve), py = 0.34 (GeV/c)^ (open squares and dotted curve), and 
Py — 0.74 (GeV/c)^ (open diamonds and dot-dashed curve). The coverage in xf is limited for the data by detector acceptance. 
The error bars are statistical only. 



hadrons [63|. In Fig. [M] the standard LO pQCD calculations are combined with a Gaussian p7-distribution (Eq. [7]), 
plotted versus xp, and compared with the data. The theory describes approximately the > behavior beginning 
from the xp peak. At negative xp values the theoretical curve is almost constant and deviates strongly from 
the data. This is because at xp < 0, z is close to zero and varies slowly, making D{z) nearly constant. In order to 
distinguish target and current fragmentation, one can use a different variable [67^ 

.c = ^ , (30) 

in which E^^'^ is hadron energy in the CM frame. This can still be interpreted as the parton momentum fraction 
carried by the measured hadron, similar to that in e^e~ collisions. By simply using the fragmentation function 
D{zg) in Eq.[2n]for both forward and backward regions, one obtains a qualitative agreement between theoretical and 
experimental xp distributions (see Fig. [25)l . Hence the target fragmentation term in Eq. [29]is equal to the standard 
"current fragmentation" contribution (1 — x)M = f{x) (g) D{zg) We speculate, therefore, that the fragmentation 
of the spectator diquark system may be quantitatively similar to the anti-quark fragmentation (see Ref. (t^ ) for 
TT"*" production. The latter mechanism is implicitly included in the fragmentation functions D{z) measured in e"'"e~ 
collisions. It is also related to the dominance of the favored u-quark fragmentation in 7r+, since the two proton's 
valence u-quarks are likely to be evenly distributed between current and target fragments. This intriguing similarity 
allows us to describe qualitatively the semi-inclusive cross section by the standard current fragmentation f(x)®D{zQ) 
term only even in the region of backward-going tt+s. 



D. Azimuthal Moments 



Figs. [26 1 [27 1 [281 and [291 show the p\ and z-dependencies oi'Hi/{'H2 + eHi) and 'Hi/{n2 + eWi). The (/)-dependent 
terms are typically less than a few percent of the ^-independent part of the semi-inclusive cross section. The (cos 20) 
moments are generally compatible with zero within our systematic uncertainties, excluding low- 2; and high-p^ region 
where they are definitely positive. The (cos (p) moments are more significant due to smaller systematic uncertainty 
and they are negative at large-pT- By exploiting the broad kinematic coverage of CLAS, we can explore the overall 
trends of the data. 

The (cos0) term shown in Fig. [21] tends to decrease as a function of pT and eventually becomes negative. For most 
of the pt range the data are 2-3 systematic deviations below zero. 
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FIG. 25: Same data as in Fig. 1241 The curves are the same as in Fig. [24] except for the fragmentation functions, which are 
evaluated at zq rather than z. 

As one can see in Fig. [27l the (cos 20) term is compatible with zero point-by-point, except for the low-z and large 
Pt where it is positive. 

The z-dependence of (cos(/)) shown in Fig. [55] has a very different behavior at the lowest px and at higher pT'- at 
the lowest px, (cos0) is compatible with zero, whereas at higher px, except for the highest px point, (cos0) rises 
from negative to positive values. 

The (cos 2(j)) term shown in Fig. [29] does not exhibit a clear z-dependence, except for low-z region where positive 
values decreasing with z can be seen. Above that region (cos 20) is generally smaller than the systematic uncertainties. 

Theoretical predictions in the (cos (f>) are in strong disagreement with our data. Indeed the full curve of the 
predictions, which has a similar dependence on px but very different z-dependence, lies many standard deviations 
below the measured points over much of the kinematics. This is due to the dominant, negative Cahn effect contribution. 
The positive contribution of the Berger effect slightly compensates for the Cahn effect, but the Berger contribution is 
too small to bring the sum of the two effects in agreement with the data. 

Theory predicts very small (cos 20) values partially due to cancellation between the Cahn and Berger effect contri- 
butions. These predictions are generally in agreement with our data. The data points at large px and low z lie above 
the theoretical curves, this difference reaches 2-3 systematic deviations. 

The averaged structure function ratios 7^3,4/(7^2 + e7ii) shown in Figs. [26l [27] [28l and [29l are listed in Tables HXl 
and |Xj We notice that the use of the weighted average technique is not strictly justified over the entire range of z 
and Px- In some points, in particular in the low-z region, the data show a clear x and/or Q^-dependence leading to 
an underestimation of the averaged statistical uncertainties. Although, the full uncertainty on the averaged data is in 
any case dominated by the systematic uncertainty. Not averaged data obtained from the two methods are statistically 
compatible. 

The comparison with higher energy data from Ref. 73] shown in Fig. [30] reveals the striking difference between the 
two measurements of H^, whereas both measurements of H4 at large and small are compatible with zero. At large 

the absolute values of the ratio H3/H2 reach 0.05-0.1 and seem to follow the expected l/Q^ behavior. However, 
our data at lower do not follow this trend having values compatible with zero. The strong suppression of H3 at 

~ 2 (GeV/c)^ with respect to the data at ~ 30 — 60 (GeV/c)^ does not seem to be related to the threshold 
effect due to the phase space shrinkage at lower energies discussed in Ref. [T^I- To account for it, the ratios of the 
Gaussian model integrals over the allowed kinematical region: 
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TABLE V: H2 structure function with its statistical and systematic uncertainties. 
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0.051 


0.50 


0.179 


0.003 


0.026 


0.07 


3.746 


0.038 


0.548 


0.41 


0.324 


0.003 


0.046 


0.45 


0.215 


0.014 


0.034 


0.54 


0.134 


0.004 


0.020 


0.09 


3.897 


0.022 


0.605 


0.45 


0.232 


0.003 


0.033 


Q^=1.74 (GoV/c)^ x=0.37 


0.59 


0.087 


0.003 


0.013 


0.12 


2.957 


0.018 


0.451 


0.50 


0.176 
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2.195 
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0.006 


0.221 
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0.026 
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0.014 
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0.004 


0.111 


0.37 


0.412 


0.007 


0.061 


0.21 


1.049 


0.005 


0.140 


0.45 


0.205 


0.010 


0.032 


0.27 


0.716 


0.003 


0.093 


0.41 


0.302 


0.007 


0.046 


0.24 


0.878 


0.005 


0.116 


0.50 


0.177 


0.009 


0.028 


0.30 


0.612 


0.003 


0.079 


0.45 


0.238 


0.007 


0.037 


0.27 


0.749 


0.005 


0.098 


Q"=1.74 (GeV/c)", x=0.41 
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0.150 


0.37 


0.389 


0.004 


0.054 


0.21 


0.878 


0.014 


0.128 


0.54 


0.122 


0.005 


0.019 


0.18 


0.980 


0.007 


0.130 


0.41 


0.302 


0.003 


0.043 


0.24 


0.726 


0.011 


0.106 


0.59 


0.097 


0.005 


0.015 


0.21 


0.897 


0.007 


0.119 


0.45 


0.220 


0.003 


0.032 


0.27 


0.608 


0.009 


0.089 


0.63 


0.071 


0.004 


0.013 


0.24 


0.702 


0.006 


0.093 


0.50 


0.176 


0.003 


0.027 


0.30 


0.535 


0.009 


0.078 


Q^=2.02 (GcV/c)^ x=0.24 


0.27 


0.595 


0.005 


0.079 


0.54 


0.139 


0.004 


0.022 


0.34 


0.434 


0.009 


0.064 


0.07 


1.683 


0.010 


0.224 


0.30 


0.524 


0.005 


0.070 


0.59 


0.127 


0.005 


0.020 


0.37 


0.356 


0.008 


0.053 


0.12 


2.002 


0.005 


0.281 


0.34 


0.412 


0.004 


0.055 


0.63 


0.083 


0.004 


0.015 


0.41 


0.308 


0.008 


0.046 


0.15 


1.617 


0.004 


0.223 


0.37 


0.352 


0.004 


0.048 


Q^=1.74 (GcV/c)^ x=0.21 


0.45 


0.258 


0.007 


0.039 


0.18 


1.326 


0.004 


0.180 


0.41 


0.314 


0.005 


0.043 


0.12 


2.106 


0.006 


0.306 


Q^=1.74 (GcV/c)^ x=0.34 


0.21 


1.115 


0.003 


0.151 


0.45 


0.244 


0.004 


0.034 


0.15 


1.717 


0.005 


0.244 


0.18 


0.867 


0.018 


0.133 


0.24 


0.915 


0.003 


0.123 


0.50 


0.213 


0.006 


0.030 


0.18 


1.383 


0.004 


0.194 


0.21 


0.788 


0.016 


0.121 


0.27 


0.774 


0.003 


0.103 


0.54 


0.165 


0.006 


0.024 


0.21 


1.155 


0.004 


0.162 


0.24 


0.640 


0.018 


0.099 


0.30 


0.632 


0.003 


0.084 


Q^=2.02 (GcV/c)^ x=0.37 


0.24 


0.949 


0.004 


0.132 


0.27 


0.569 


0.012 


0.088 


0.34 


0.514 


0.003 


0.069 


0.15 


0.950 


0.010 


0.128 


0.27 


0.788 


0.003 


0.109 


0.30 


0.487 


0.015 


0.076 


0.37 


0.410 


0.003 


0.055 


0.18 


0.861 


0.008 


0.115 


0.30 


0.644 


0.003 


0.089 


0.34 


0.363 


0.010 


0.056 


0.41 


0.312 


0.003 


0.043 


0.21 


0.784 


0.008 


0.105 


0.34 


0.518 


0.003 


0.072 


0.37 


0.330 


0.014 


0.053 


0.45 


0.237 


0.003 


0.033 


0.24 


0.608 


0.007 


0.082 


0.27 


0.528 


0.005 


0.072 


0.54 


0.105 


0.006 


0.016 


0.34 


0.384 


0.004 


0.049 


0.15 


1.645 


0.008 


0.230 


0.30 


0.488 


0.006 


0.066 


Q^=2.37 (GeV/c)^ x=0.27 


0.37 


0.308 


0.004 


0.040 


0.18 


1.345 


0.006 


0.186 



TABLE VI: Continue from Table |Vl 



2 


H2 


stat.err. 


sys.err. 


z 


H2 


stat.err. 


sys.err. 


z 


Hi 


stat.err. 


sys.err. 


2 


^2 


stat.err. 


sys.err. 


0.34 


0.369 


0.005 


0.050 


0.07 


2.047 


0.014 


0.280 


0.41 


0.286 


0.004 


0.037 


0.21 


1.113 


0.006 


0.153 


0.37 


0.320 


0.005 


0.044 


0.12 


2.162 


0.006 


0.311 


0.45 


0.228 


0.004 


0.030 


0.24 


0.926 


0.006 


0.127 


0.41 


0.300 


0.006 


0.042 


0.15 


1.730 


0.005 


0.244 


0.50 


0.187 


0.005 


0.025 


0.27 


0.761 


0.005 


0.104 


0.45 


0.214 


0.005 


0.030 


0.18 


1.437 


0.005 


0.200 


0.54 


0.153 


0.005 


0.021 


0.30 


0.621 


0.005 


0.085 


0.50 


0.188 


0.005 


0.027 


0.21 


1.196 


0.004 


0.166 


Q^=2.37 (GeV/c)", x=0.41 


0.34 


0.505 


0.006 


0.069 


Q^=2.02 (GeV/c)^ x=0.41 


0.24 


0.998 


0.004 


0.138 


0.15 


0.895 


0.008 


0.113 


0.37 


0.388 


0.005 


0.054 


0.18 


0.700 


0.012 


0.102 


0.27 


0.824 


0.004 


0.113 


0.18 


0.789 


0.006 


0.099 


0.41 


0.339 


0.007 


0.048 


0.21 


0.639 


0.008 


0.094 


0.30 


0.675 


0.004 


0.093 


0.21 


0.663 


0.007 


0.084 


0.45 


0.224 


0.005 


0.032 


0.24 


0.495 


0.009 


0.073 


0.34 


0.562 


0.004 


0.077 


0.24 


0.580 


0.005 


0.073 


Q^=2.93 (GeV/c)^ x=0.34 


0.27 


0.449 


0.007 


0.066 


0.37 


0.440 


0.004 


0.061 


0.27 


0.477 


0.004 


0.061 


0.07 


1.469 


0.016 


0.191 


0.30 


0.420 


0.007 


0.063 


0.41 


0.341 


0.004 


0.048 


0.30 


0.404 


0.004 


0.052 


0.09 


1.926 


0.007 


0.268 


0.34 


0.359 


0.006 


0.054 


0.45 


0.251 


0.004 


0.036 


0.34 


0.347 


0.005 


0.045 


0.12 


1.629 


0.006 


0.225 


0.37 


0.293 


0.005 


0.045 


0.50 


0.188 


0.004 


0.028 


0.37 


0.277 


0.004 


0.036 


0.15 


1.378 


0.005 


0.187 


0.41 


0.256 


0.005 


0.039 


0.54 


0.137 


0.004 


0.021 


0.41 


0.257 


0.005 


0.034 


0.18 


1.124 


0.004 


0.150 


0.45 


0.216 


0.004 


0.033 


0.59 


0.101 


0.004 


0.016 


0.45 


0.200 


0.005 


0.027 


0.21 


0.943 


0.004 


0.125 


0.50 


0.188 


0.004 


0.029 


0.68 


0.081 


0.020 


0.016 


0.50 


0.153 


0.003 


0.021 


0.24 


0.804 


0.004 


0.106 


0.54 


0.147 


0.004 


0.023 


Q^=2.37 (GeV/c)^ x=0.30 


Q"=2.37 (GeV/c)", x=0.45 


0.27 


0.672 


0.004 


0.088 


Q^=2.02 (GeV/c)^ x=0.45 


0.15 


1.401 


0.005 


0.185 


0.15 


0.898 


0.012 


0.120 


0.30 


0.560 


0.004 


0.073 


0.18 


0.721 


0.011 


0.116 


0.18 


1.134 


0.004 


0.147 


0.18 


0.628 


0.009 


0.084 


0.34 


0.484 


0.004 


0.063 


0.21 


0.489 


0.009 


0.079 


0.21 


0.989 


0.004 


0.128 


0.21 


0.553 


0.006 


0.074 


0.37 


0.370 


0.004 


0.049 


0.24 


0.392 


0.006 


0.063 


0.24 


0.852 


0.004 


0.109 


0.24 


0.469 


0.006 


0.063 


0.41 


0.283 


0.005 


0.038 


0.27 


0.368 


0.007 


0.060 


0.27 


0.704 


0.004 


0.090 


0.27 


0.410 


0.006 


0.056 


0.45 


0.221 


0.004 


0.030 


0.30 


0.348 


0.006 


0.057 


0.30 


0.590 


0.003 


0.075 


0.30 


0.357 


0.005 


0.049 


0.50 


0.145 


0.005 


0.021 


0.34 


0.295 


0.006 


0.049 


0.34 


0.511 


0.004 


0.065 


0.34 


0.309 


0.005 


0.043 


0.54 


0.110 


0.004 


0.016 


0.37 


0.247 


0.005 


0.041 


0.37 


0.395 


0.004 


0.050 


0.37 


0.246 


0.004 


0.035 


Q^=2.93 (GeV/c)^ x=0.37 


0.41 


0.226 


0.005 


0.037 


0.41 


0.316 


0.004 


0.041 


0.41 


0.234 


0.004 


0.033 


0.09 


1.333 


0.013 


0.174 


0.45 


0.188 


0.006 


0.031 


0.45 


0.233 


0.004 


0.030 


0.50 


0.168 


0.003 


0.024 


0.12 


1.282 


0.006 


0.168 


Q^=2.02 (GeV/c)^ x=0.50 


0.50 


0.156 


0.004 


0.021 


Q"=2.37 (GeV/c)", x=0.50 


0.15 


1.138 


0.006 


0.149 


0.21 


0.519 


0.011 


0.076 


0.54 


0.152 


0.004 


0.021 


0.15 


0.492 


0.034 


0.065 


0.18 


0.948 


0.005 


0.122 


0.24 


0.357 


0.008 


0.053 


0.59 


0.115 


0.005 


0.016 


0.18 


0.590 


0.007 


0.078 


0.21 


0.785 


0.004 


0.100 


0.27 


0.289 


0.010 


0.043 


0.63 


0.095 


0.006 


0.014 


0.21 


0.435 


0.007 


0.059 


0.24 


0.711 


0.004 


0.090 


0.30 


0.300 


0.008 


0.045 


Q^=2.37 (GeV/c)^ x=0.34 


0.24 


0.394 


0.005 


0.053 


0.27 


0.609 


0.004 


0.077 


0.34 


0.246 


0.006 


0.037 


0.09 


1.094 


0.013 


0.132 


0.27 


0.350 


0.005 


0.048 


0.30 


0.505 


0.005 


0.063 


0.37 


0.197 


0.006 


0.030 


0.12 


1.295 


0.006 


0.157 


0.30 


0.305 


0.005 


0.042 


0.34 


0.427 


0.004 


0.054 


0.41 


0.181 


0.007 


0.027 


0.15 


1.197 


0.006 


0.146 


0.34 


0.295 


0.005 


0.041 


0.37 


0.316 


0.004 


0.040 


Q^=2.02 (GeV/c)^ x=0.54 


0.18 


0.990 


0.005 


0.120 


0.37 


0.226 


0.004 


0.032 


0.41 


0.250 


0.004 


0.032 


0.18 


0.162 


0.015 


0.026 


0.21 


0.832 


0.004 


0.101 


Q"=2.37 (GeV/c)", x=0.54 


0.45 


0.253 


0.005 


0.033 


0.21 


0.207 


0.007 


0.034 


0.24 


0.742 


0.004 


0.090 


0.18 


0.316 


0.024 


0.042 


0.50 


0.171 


0.004 


0.023 


0.24 


0.175 


0.009 


0.029 


0.27 


0.622 


0.004 


0.076 


0.21 


0.369 


0.006 


0.049 


0.54 


0.172 


0.005 


0.024 


0.27 


0.176 


0.006 


0.029 


0.30 


0.544 


0.005 


0.067 


0.24 


0.332 


0.014 


0.045 


0.59 


0.117 


0.005 


0.016 


Q^=2.37 (GeV/c)^ x=0.24 


0.34 


0.437 


0.004 


0.054 


0.27 


0.297 


0.009 


0.041 


Q^=2.93 (GeV/c)^ x=0.41 


0.07 


3.085 


0.030 


0.419 


0.37 


0.343 


0.004 


0.043 


0.30 


0.246 


0.005 


0.033 


0.12 


1.089 


0.007 


0.128 


0.12 


2.502 


0.014 


0.354 


0.41 


0.288 


0.004 


0.036 


0.34 


0.229 


0.005 


0.031 


0.15 


0.961 


0.006 


0.114 


0.15 


1.860 


0.012 


0.258 


0.45 


0.220 


0.003 


0.028 


0.37 


0.185 


0.005 


0.026 


0.18 


0.795 


0.007 


0.094 


0.18 


1.477 


0.010 


0.202 


0.50 


0.182 


0.004 


0.024 


Q"=2.37 (GeV/c)", x=0.59 


0.21 


0.668 


0.005 


0.079 


0.21 


1.224 


0.009 


0.168 


0.54 


0.143 


0.003 


0.019 


0.18 


0.175 


0.007 


0.024 


0.24 


0.600 


0.005 


0.072 


0.24 


1.000 


0.008 


0.136 


Q^=2.37 (GeV/c)^ x=0.37 


0.21 


0.180 


0.004 


0.025 


0.27 


0.544 


0.005 


0.065 


0.27 


0.805 


0.008 


0.109 


0.09 


0.755 


0.026 


0.093 


0.24 


0.162 


0.004 


0.023 


0.30 


0.418 


0.004 


0.051 


0.30 


0.653 


0.008 


0.089 


0.15 


1.044 


0.008 


0.130 


0.27 


0.142 


0.005 


0.020 


0.34 


0.367 


0.004 


0.045 


0.34 


0.509 


0.008 


0.070 


0.18 


0.886 


0.007 


0.110 


0.30 


0.136 


0.005 


0.019 


0.37 


0.277 


0.004 


0.034 


0.37 


0.387 


0.007 


0.054 


0.21 


0.762 


0.006 


0.095 


0.34 


0.129 


0.003 


0.018 


0.41 


0.235 


0.004 


0.029 


0.41 


0.306 


0.006 


0.043 


0.24 


0.653 


0.005 


0.082 


Q"=2.93 (GeV/c)", x=0.30 


0.45 


0.240 


0.006 


0.031 


0.45 


0.221 


0.006 


0.032 


0.27 


0.540 


0.004 


0.068 


0.09 


2.623 


0.010 


0.377 


0.50 


0.167 


0.005 


0.022 


0.50 


0.160 


0.006 


0.024 


0.30 


0.477 


0.005 


0.061 


0.12 


2.089 


0.008 


0.297 


Q^=2.93 (GeV/c)^ x=0.45 


0.09 


0.522 


0.024 


0.063 


0.24 


0.840 


0.018 


0.117 


0.24 


0.435 


0.005 


0.054 


0.24 


0.566 


0.008 


0.076 


0.12 


0.940 


0.010 


0.113 


0.27 


0.655 


0.016 


0.091 


0.27 


0.383 


0.004 


0.048 


0.27 


0.457 


0.008 


0.061 


0.15 


0.794 


0.007 


0.096 


0.30 


0.545 


0.017 


0.076 


0.30 


0.313 


0.005 


0.039 


0.30 


0.416 


0.009 


0.056 


0.18 


0.738 


0.007 


0.089 


0.34 


0.445 


0.017 


0.062 


0.34 


0.259 


0.004 


0.032 


0.34 


0.329 


0.009 


0.044 



TABLE VIL Continue from Table El 



2 


H2 


stat.err. 


sys.err. 


z 




stat.err. 


sys.err. 


z 


Hi 


stat.err. 


sys.err. 


2 


^2 


stat.err. 


sys.err. 


0.21 


0.560 


0.007 


0.068 


0.37 


0.335 


0.015 


0.047 


0.37 


0.247 


0.006 


0.032 


0.37 


0.231 


0.011 


0.031 


0.24 


0.509 


0.005 


0.062 


0.41 


0.328 


0.022 


0.047 


0.41 


0.184 


0.005 


0.024 


0.41 


0.206 


0.011 


0.028 


0.27 


0.461 


0.005 


0.057 


0.45 


0.199 


0.020 


0.029 


0.45 


0.182 


0.005 


0.024 


0.45 


0.138 


0.013 


0.019 


0.30 


0.342 


0.004 


0.042 


Q^=3.42 (GeV/c)^ x=0.37 


0.50 


0.132 


0.004 


0.018 


Q^=4.10 (GeV/c)^ x=0.50 


0.34 


0.300 


0.004 


0.037 


0.09 


1.847 


0.011 


0.261 


Q"=3.42 (GeV/c)", x=0.54 


0.09 


0.950 


0.010 


0.134 


0.37 


0.236 


0.005 


0.030 


0.12 


1.521 


0.008 


0.214 


0.12 


0.775 


0.012 


0.091 


0.12 


0.872 


0.007 


0.121 


0.41 


0.215 


0.005 


0.027 


0.15 


1.242 


0.008 


0.172 


0.15 


0.516 


0.007 


0.061 


0.15 


0.747 


0.008 


0.102 


0.45 


0.213 


0.005 


0.028 


0.18 


1.033 


0.006 


0.141 


0.18 


0.466 


0.006 


0.055 


0.18 


0.609 


0.006 


0.082 


0.50 


0.133 


0.003 


0.018 


0.21 


0.880 


0.006 


0.119 


0.21 


0.355 


0.006 


0.043 


0.21 


0.546 


0.007 


0.073 


0.54 


0.143 


0.003 


0.020 


0.24 


0.749 


0.006 


0.101 


0.24 


0.321 


0.004 


0.039 


0.24 


0.481 


0.006 


0.064 


Q^=2.93 (GeV/c)^ x=0.50 


0.27 


0.617 


0.006 


0.083 


0.27 


0.335 


0.005 


0.041 


0.27 


0.393 


0.006 


0.052 


0.15 


0.734 


0.010 


0.089 


0.30 


0.525 


0.006 


0.070 


0.30 


0.247 


0.004 


0.031 


0.30 


0.357 


0.007 


0.047 


0.18 


0.619 


0.007 


0.075 


0.34 


0.442 


0.006 


0.059 


0.34 


0.206 


0.004 


0.026 


0.34 


0.276 


0.006 


0.037 


0.21 


0.432 


0.007 


0.053 


0.37 


0.349 


0.006 


0.047 


0.37 


0.221 


0.004 


0.028 


0.37 


0.214 


0.006 


0.029 


0.24 


0.405 


0.006 


0.050 


0.41 


0.258 


0.007 


0.035 


0.41 


0.159 


0.003 


0.020 


0.41 


0.180 


0.006 


0.024 


0.27 


0.457 


0.007 


0.058 


0.45 


0.198 


0.007 


0.028 


0.45 


0.151 


0.004 


0.020 


0.45 


0.135 


0.007 


0.019 


0.30 


0.306 


0.005 


0.039 


Q^=3.42 (GeV/c)^ x=0.41 


Q"=3.42 (GeV/c)", x=0.59 


Q^=4.10 (GeV/c)^ x=0.54 


0.34 


0.281 


0.005 


0.036 


0.09 


1.458 


0.009 


0.200 


0.15 


0.406 


0.006 


0.049 


0.09 


0.539 


0.017 


0.066 


0.37 


0.219 


0.003 


0.028 


0.12 


1.270 


0.007 


0.173 


0.18 


0.356 


0.006 


0.043 


0.12 


0.629 


0.007 


0.077 


0.41 


0.195 


0.003 


0.025 


0.15 


1.025 


0.007 


0.137 


0.21 


0.291 


0.007 


0.036 


0.15 


0.571 


0.007 


0.070 


0.50 


0.135 


0.003 


0.018 


0.18 


0.866 


0.005 


0.114 


0.24 


0.259 


0.004 


0.032 


0.18 


0.458 


0.006 


0.056 


Q^=2.93 (GeV/c)^ x=0.54 


0.21 


0.770 


0.005 


0.101 


0.27 


0.237 


0.004 


0.030 


0.21 


0.426 


0.006 


0.052 


0.15 


0.605 


0.008 


0.074 


0.24 


0.662 


0.005 


0.086 


0.30 


0.204 


0.004 


0.026 


0.24 


0.364 


0.005 


0.045 


0.18 


0.491 


0.007 


0.060 


0.27 


0.548 


0.005 


0.071 


0.34 


0.167 


0.003 


0.021 


0.27 


0.292 


0.006 


0.036 


0.21 


0.327 


0.006 


0.040 


0.30 


0.474 


0.006 


0.061 


Q^=3.42 (GeV/c)", x=0.63 


0.30 


0.268 


0.006 


0.033 


0.24 


0.322 


0.005 


0.040 


0.34 


0.382 


0.005 


0.049 


0.15 


0.193 


0.007 


0.024 


0.34 


0.233 


0.005 


0.029 


0.27 


0.329 


0.005 


0.042 


0.37 


0.300 


0.005 


0.039 


0.18 


0.233 


0.004 


0.029 


0.37 


0.204 


0.007 


0.026 


0.30 


0.242 


0.004 


0.031 


0.41 


0.231 


0.005 


0.030 


0.21 


0.159 


0.006 


0.020 


0.41 


0.162 


0.006 


0.021 


0.34 


0.221 


0.004 


0.028 


0.45 


0.200 


0.006 


0.026 


0.24 


0.187 


0.009 


0.024 


0.45 


0.130 


0.006 


0.017 


0.37 


0.203 


0.003 


0.026 


0.50 


0.163 


0.006 


0.022 


0.27 


0.181 


0.007 


0.023 


0.50 


0.110 


0.004 


0.015 


0.41 


0.136 


0.003 


0.018 


0.54 


0.128 


0.006 


0.018 


0.30 


0.155 


0.006 


0.020 


Q^=4.10 (GeV/c)^ x=0.59 


Q^=2.93 (GeV/c)^ x=0.59 


Q^=3.42 (GeV/c)^ x=0.45 


0.34 


0.117 


0.004 


0.015 


0.09 


0.408 


0.015 


0.047 


0.18 


0.396 


0.006 


0.052 


0.09 


0.982 


0.009 


0.118 


Q"=3.42 (GeV/c)", x=0.68 


0.12 


0.510 


0.007 


0.059 


0.21 


0.280 


0.005 


0.037 


0.12 


0.978 


0.007 


0.117 


0.15 


0.083 


0.004 


0.010 


0.15 


0.427 


0.008 


0.050 


0.24 


0.248 


0.007 


0.033 


0.15 


0.848 


0.007 


0.102 


0.18 


0.104 


0.003 


0.012 


0.18 


0.337 


0.006 


0.039 


0.27 


0.224 


0.007 


0.030 


0.18 


0.735 


0.006 


0.089 


0.24 


0.082 


0.003 


0.010 


0.21 


0.298 


0.006 


0.036 


0.30 


0.177 


0.004 


0.024 


0.21 


0.623 


0.005 


0.075 


0.27 


0.088 


0.003 


0.011 


0.24 


0.312 


0.007 


0.037 


0.34 


0.162 


0.005 


0.022 


0.24 


0.534 


0.005 


0.065 


Q"=4.10 (GeV/c)", x=0.41 


0.27 


0.249 


0.005 


0.030 


0.37 


0.149 


0.005 


0.020 


0.27 


0.478 


0.005 


0.058 


0.09 


1.718 


0.049 


0.259 


0.30 


0.229 


0.005 


0.028 


Q^=2.93 (GeV/c)^ x=0.63 


0.30 


0.382 


0.005 


0.047 


0.12 


1.406 


0.039 


0.209 


0.34 


0.176 


0.004 


0.022 


0.15 


0.122 


0.006 


0.017 


0.34 


0.308 


0.004 


0.038 


0.15 


0.992 


0.034 


0.146 


0.37 


0.144 


0.004 


0.018 


0.18 


0.154 


0.004 


0.022 


0.37 


0.263 


0.005 


0.033 


0.18 


0.918 


0.029 


0.133 


0.41 


0.154 


0.004 


0.019 


0.21 


0.138 


0.003 


0.020 


0.41 


0.211 


0.004 


0.026 


0.21 


0.764 


0.030 


0.110 


Q^=4.10 (GeV/c)^ x=0.63 


0.24 


0.125 


0.003 


0.018 


0.45 


0.212 


0.007 


0.028 


0.24 


0.628 


0.026 


0.090 


0.15 


0.333 


0.006 


0.041 


0.27 


0.150 


0.005 


0.022 


0.50 


0.150 


0.006 


0.020 


0.30 


0.383 


0.024 


0.055 


0.18 


0.264 


0.006 


0.032 


0.34 


0.099 


0.002 


0.014 


0.54 


0.137 


0.007 


0.019 


Q^=4.10 (GeV/c)", x=0.45 


0.21 


0.218 


0.005 


0.027 


Q^=3.42 (GeV/c)^ x=0.34 


Q^=3.42 (GeV/c)^ x=0.50 


0.07 


1.033 


0.033 


0.137 


0.24 


0.207 


0.004 


0.026 


0.07 


2.260 


0.064 


0.312 


0.09 


0.542 


0.008 


0.066 


0.09 


1.245 


0.013 


0.177 


0.27 


0.168 


0.004 


0.022 


0.12 


1.835 


0.024 


0.264 


0.12 


0.827 


0.010 


0.101 


0.12 


1.102 


0.011 


0.156 


0.30 


0.168 


0.004 


0.022 


0.15 


1.379 


0.021 


0.196 


0.15 


0.670 


0.006 


0.082 


0.15 


0.896 


0.010 


0.124 


0.34 


0.141 


0.004 


0.018 


0.18 


1.133 


0.018 


0.159 


0.18 


0.592 


0.007 


0.072 


0.18 


0.746 


0.008 


0.102 


Q^=4.10 (GeV/c)^ x=0.68 


0.21 


0.977 


0.018 


0.137 


0.21 


0.510 


0.006 


0.063 


0.21 


0.648 


0.008 


0.088 


0.15 


0.147 


0.003 


0.018 


0.18 


0.131 


0.003 


0.016 


0.09 


0.599 


0.012 


0.081 


0.15 


0.318 


0.008 


0.042 


0.30 


0.108 


0.004 


0.014 


0.21 


0.136 


0.006 


0.017 


0.12 


0.492 


0.010 


0.064 


0.18 


0.228 


0.007 


0.030 


0.34 


0.100 


0.004 


0.013 


0.24 


0.131 


0.006 


0.016 


0.15 


0.435 


0.010 


0.057 


0.21 


0.216 


0.006 


0.029 


Q^=4.85 (GeV/c)^ x=0.74 


Q^=4.10 (GeV/c)^ x=0.74 


0.18 


0.370 


0.008 


0.048 


0.24 


0.189 


0.006 


0.025 


0.15 


0.086 


0.003 


0.010 


0.15 


0.076 


0.005 


0.009 


0.21 


0.321 


0.008 


0.042 


0.27 


0.170 


0.005 


0.023 


0.21 


0.082 


0.004 


0.010 


0.18 


0.067 


0.003 


0.008 


0.24 


0.304 


0.008 


0.040 


0.30 


0.138 


0.006 


0.019 


Q^=5.72 (GeV/c)^ x=0.68 
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TABLE VIII: Continue from Table El 



2 1 Hi 


stat.err. | sys.err. 


z 


Hi 


stat.err. 


sys.err. 


z 


Hi 


stat.err. 


sys.err. 


2 


Hi 


stat.err. 


sys.err. 


Q^=4.85 (GeV/c)^ x=0.54 


0.27 


0.230 


0.007 


0.030 


0.34 


0.144 


0.005 


0.020 


0.09 


0.308 


0.024 


0.047 


0.09 


0.817 


0.023 


0.124 


0.30 


0.202 


0.009 


0.027 


0.37 


0.111 


0.005 


0.016 


0.12 


0.265 


0.011 


0.041 


0.12 


0.658 


0.016 


0.098 


0.34 


0.184 


0.011 


0.025 


Q"=4.85 (GeV/c)", x=0.68 


0.15 


0.205 


0.014 


0.032 


0.15 


0.600 


0.018 


0.088 


0.37 


0.127 


0.009 


0.017 


0.15 


0.218 


0.005 


0.027 


0.18 


0.156 


0.013 


0.024 


0.18 


0.499 


0.014 


0.073 


0.41 


0.141 


0.013 


0.019 


0.18 


0.227 


0.009 


0.028 


0.21 


0.124 


0.010 


0.019 


0.21 


0.424 


0.014 


0.061 


Q^=4.85 (GeV/c)^ x=0.63 


0.21 


0.145 


0.006 


0.018 


0.30 


0.114 


0.011 


0.018 


0.27 


0.305 


0.015 


0.044 


0.09 


0.345 


0.012 


0.045 


0.24 


0.136 


0.004 


0.017 


Q^=5.72 (GeV/c)^ x=0.74 


0^=4.85 (GeV/c)^ x=0.59 


0.12 


0.350 


0.006 


0.046 


0.27 


0.129 


0.004 


0.017 


0.15|0.127| 0.008 1 0.019 



for H'i/H2 and 

/(^"■l^ pleyi^[-pl/{pl)]dpl 

for H4/H2 are included in the Cahn effect curves in Fig. [3D1 These corrections do not affect strongly the Cahn effect 
curves in the presented interval of Q^. 



V. CONCLUSIONS 



We performed a measurement of semi-inclusive tt^ electroproduction in the range from 1.4 to 5.7 (GeV/c)^ 
with broad coverage in all other kinematic variables. The five-fold differential cross sections allowed us to separate 
the contributions of different structure functions. From these data we draw the following conclusions: 

• The transverse momentum dependence for the ^-independent term Ti.2 + tHi exhibits the expected thermal 
Gaussian distribution. 

• At large z the mean transverse momentum {p^) is found to be x and Q^-independent and it rises with z as 
expected within the naive parton model. In the low- 2: region, {p^) is altered by the limited phase space. 

• The comparison of the measured structure function H2 to the current fragmentation LO and NLO pQCD 
calculations shows that the difference between the data and calculations reaches 20% at the low- a; limit of 
the z < 0.4 region, which is compatible with systematic uncertainties in the calculations due to higher-order 
corrections and the favored fragmentation assumption. 

• The separation of the current and soft target fragmentation processes in the CM frame shows a symmetry 
about xf = 0, which suggests the presence of an intriguing numerical equality between the fragmentation of the 
spectator diquark in the target region of SIDIS and the fragmentation of the anti-quark in e+e~ collisions. 

• The precision of the data does not allow us to obtain information about the contribution of the Boer-Mulders 
function, which is expected to be smaller than our estimated systematic uncertainties js^ . [75| . 

• The 7^4 structure function is compatible with zero within our precision, except for the low- 2; region where it 
is positive. The Ti.^ structure function appears to be somewhat better determined than Ti.^ and is in strong 
disagreement with the predictions of the Cahn effect. Inclusion of the Berger effect does not change significantly 
the disagreement in H3. 

• H3 structure function at ~ 2 (GcV/c)^ found to be strongly suppressed with respect to the data at ~ 
30 — 60 (GeV/c)^. This suppression does not seem to be related to the the phase space shrinkage at our energies. 



The data tables can be found in the CLAS physics database [t^ and in sources of the electronic preprint [tJ 
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FIG. 26: The p^-dependence of 7^3/(7^2 + eTYi) (open squares - moments, full triangles - fits) for different z averaged over x 
and . The thick curves show theoretical predictions of the Cahn effect [3ll. Is^ (dashed), predictions of the Berger effect [4l[ 
using a convex pion wave function (dotted) and their sum (solid). The two data sets (from moments and fits extractions) 
are shifted equally along the x-axis in opposite directions from their central values for visibility. The hatched area show the 
systematic uncertainties. 
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FIG. 27: Same as Fig. 1261 except for 7^4/(7^2 + eTii) arising from the {cos2(j)) moment (open triangles - moments, full squares 
- fits). 
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TABLE IX: The extracted data on 7^3/(7-^2 + eHi) averaged over x and with their statistical and systematic uncertainties. 
The values given in brackets are obtained using the fit method. 



z 


Pt 
{GcV/cf 


(v) 


{GeV/cf 


(x) 


Hs/iHi+eHi) 


stat. uncertainty 


sys. uncertainty 


0.068 


0.026 


0.88 (0.87) 


2.12 (2.22) 


0.23 (0.24) 


0.012 ( 0.012) 


0.0007 (0.0007) 


0.005 


0.068 


0.065 


0.88 (0.87) 


2.14 (2.19) 


0.23 (0.23) 


0.001 (-0.001) 


0.0004 (0.0005) 


0.008 


0.093 


0.026 


0.85 (0.84) 


2.14 (2.30) 


0.24 (0.26) 


0.009 ( 0.010) 


0.0005 (0.0005) 


0.006 


0.093 


0.065 


0.83 (0.83) 


2.27 (2.31) 


0.26 (0.26) 


-0.001 (-0.003) 


0.0002 (0.0003) 


0.007 


0.093 


0.129 


0.87 (0.87) 


2.10 (2.13) 


0.22 (0.23) 


-0.004 (-0.008) 


0.0004 (0.0004) 


0.014 


0.093 


0.217 


0.92 (0.92) 


1.85 (1.86) 


0.19 (0.19) 


-0.021 (-0.023) 


0.0015 (0.0017) 


0.002 


0.119 


0.026 


0.84 (0.83) 


2.12 (2.26) 


0.24 (0.26) 


0.006 ( 0.009) 


0.0009 (0.0008) 


0.005 


0.119 


0.065 


0.81 (0.81) 


2.28 (2.35) 


0.27 (0.27) 


-0.004 (-0.005) 


0.0003 (0.0003) 


0.008 


0.119 


0.129 


0.84 (0.84) 


2.23 (2.27) 


0.25 (0.25) 


-0.003 (-0.008) 


0.0003 (0.0003) 


0.006 


0.119 


0.217 


0.88 (0.88) 


2.01 (2.03) 


0.21 (0.22) 


-0.010 (-0.013) 


0.0006 (0.0006) 


0.007 


0.119 


0.340 


0.93 (0.93) 


1.66 (1.66) 


0.17 (0.17) 


-0.025 (-0.027) 


0.0033 (0.0035) 


0.002 


0.147 


0.005 


0.61 (0.59) 


2.83 (2.65) 


0.42 (0.41) 


0.020 (-0.040) 


0.0298 (0.0107) 


0.013 


0.147 


0.026 


0.89 (0.91) 


1.89 (1.99) 


0.20 (0.20) 


0.007 ( 0.012) 


0.0035 (0.0013) 


0.003 


0.147 


0.065 


0.81 (0.80) 


2.20 (2.32) 


0.26 (0.27) 


-0.003 (-0.003) 


0.0004 (0.0004) 


0.008 


0.147 


0.129 


0.82 (0.82) 


2.23 (2.28) 


0.26 (0.26) 


-0.009 (-0.011) 


0.0003 (0.0004) 


0.007 


0.147 


0.217 


0.85 (0.85) 


2.14 (2.18) 


0.24 (0.24) 


-0.008 (-0.013) 


0.0004 (0.0005) 


0.005 


0.147 


0.340 


0.91 (0.91) 


1.99 (2.01) 


0.20 (0.20) 


-0.014 (-0.016) 


0.0011 (0.0013) 


0.005 


0.175 


0.005 


0.59 (0.60) 


2.57 (2.55) 


0.40 (0.39) 


0.009 ( 0.004) 


0.0109 (0.0052) 


0.006 


0.175 


0.065 


0.82 (0.81) 


2.19 (2.31) 


0.25 (0.27) 


0.000 ( 0.000) 


0.0005 (0.0005) 


0.004 


0.175 


0.129 


0.82 (0.82) 


2.28 (2.34) 


0.26 (0.27) 


-0.010 (-0.011) 


0.0004 (0.0004) 


0.006 


0.175 


0.217 


0.84 (0.84) 


2.20 (2.24) 


0.24 (0.25) 


-0.011 (-0.013) 


0.0004 (0.0005) 


0.005 


0.175 


0.340 


0.87 (0.87) 


2.02 (2.04) 


0.22 (0.22) 


-0.012 (-0.016) 


0.0007 (0.0008) 


0.006 


0.175 


0.507 


0.92 (0.92) 


1.85 (1.86) 


0.19 (0.19) 


-0.010 (-0.011) 


0.0024 (0.0026) 


0.008 


0.205 


0.005 


0.60 (0.62) 


2.65 (2.68) 


0.41 (0.40) 


0.009 ( 0.008) 


0.0056 (0.0035) 


0.006 


0.205 


0.065 


0.82 (0.81) 


2.12 (2.23) 


0.25 (0.26) 


0.010 ( 0.008) 


0.0008 (0.0007) 


0.004 


0.205 


0.129 


0.81 (0.81) 


2.27 (2.34) 


0.26 (0.27) 


-0.006 (-0.008) 


0.0004 (0.0005) 


0.005 


0.205 


0.217 


0.84 (0.83) 


2.24 (2.28) 


0.25 (0.26) 


-0.012 (-0.013) 


0.0005 (0.0005) 


0.005 


0.205 


0.340 


0.87 (0.86) 


2.14 (2.17) 


0.23 (0.23) 


-0.011 (-0.012) 


0.0007 (0.0008) 


0.006 


0.205 


0.507 


0.92 (0.92) 


1.85 (1.87) 


0.19 (0.19) 


-0.006 (-0.010) 


0.0018 (0.0020) 


0.007 


0.236 


0.005 


0.59 (0.61) 


2.56 (2.61) 


0.40 (0.39) 


0.013 ( 0.008) 


0.0036 (0.0027) 


0.006 


0.236 


0.065 


0.86 (0.86) 


1.86 (1.87) 


0.20 (0.20) 


0.013 ( 0.013) 


0.0024 (0.0020) 


0.004 


0.236 


0.129 


0.81 (0.81) 


2.23 (2.32) 


0.26 (0.27) 


0.001 (-0.002) 


0.0005 (0.0005) 


0.006 


0.236 


0.217 


0.82 (0.82) 


2.23 (2.28) 


0.26 (0.26) 


-0.011 (-0.013) 


0.0005 (0.0005) 


0.005 


0.236 


0.340 


0.84 (0.84) 


2.15 (2.17) 


0.24 (0.24) 


-0.016 (-0.017) 


0.0006 (0.0007) 


0.006 


0.236 


0.507 


0.90 (0.90) 


1.97 (1.99) 


0.20 (0.21) 


-0.009 (-0.010) 


0.0014 (0.0016) 


0.010 


0.236 


0.7 l:i 


().9;i (().9;i) 


i.()() (i.()7) 


0.17 (0.17) 


0.001 (-0.002) 


O.OOoO (O.OO.IS) 


0.001 


0.269 


0.005 


0.61 (0.63) 


2.66 (2.67) 


0.40 (0.39) 


0.011 ( 0.007) 


0.0033 (0.0027) 


0.006 


0.269 


0.129 


0.81 (0.83) 


2.17 (2.14) 


0.25 (0.24) 


0.009 ( 0.005) 


0.0007 (0.0007) 


0.006 


0.269 


0.217 


0.81 (0.81) 


2.23 (2.28) 


0.26 (0.26) 


-0.008 (-0.010) 


0.0005 (0.0006) 


0.005 


0.269 


o.;mo 


0.84 (0.84) 


2.15 (2.19) 


0.24 (0.24) 


-0.018 (-0.019) 


0.0007 (0.0008) 


0.007 


0.269 


0..j()7 


0.87 (0.87) 


2.():i (2.()()) 


0.22 (0.22) 


-0.011 (-0.015) 


O.OOll (0.001:5) 


0.009 


0.269 


0.743 


0.92 (0.92) 


1.77 (1.78) 


0.18 (0.18) 


0.005 ( 0.003) 


0.0036 (0.0041) 


0.008 


0.303 


0.005 


0.62 (0.64) 


2.48 (2.53) 


0.37 (0.37) 


0.012 ( 0.014) 


0.0035 (0.0029) 


0.007 


0.303 


0.129 


0.83 (0.82) 


2.03 (2.11) 


0.23 (0.24) 


0.012 ( 0.010) 


0.0012 (0.0012) 


0.006 


0.303 


0.217 


0.82 (0.82) 


2.21 (2.29) 


0.25 (0.26) 


-0.002 (-0.004) 


0.0006 (0.0007) 


0.005 


0.303 


0.340 


0.83 (0.84) 


2.16 (2.20) 


0.24 (0.25) 


-0.015 (-0.016) 


0.0007 (0.0008) 


0.005 


0.303 


0.507 


0.87 (0.87) 


2.03 (2.05) 


0.22 (0.22) 


-0.021 (-0.022) 


0.0011 (0.0013) 


0.008 


0.303 


0.743 


0.92 (0.92) 


1.86 (1.86) 


0.19 (0.19) 


0.003 ( 0.000) 


0.0029 (0.0033) 


0.009 


0.338 


0.005 


0.63 (0.64) 


2.50 (2.53) 


0.37 (0.36) 


0.016 ( 0.015) 


0.0033 (0.0029) 


0.007 


0.338 


0.026 


0.63 (0.62) 


2.71 (2.68) 


0.40 (0.40) 


0.002 (-0.009) 


0.0386 (0.0138) 


0.006 


0.338 


0.129 


0.92 (0.92) 


1.74 (1.74) 


0.18 (0.18) 


0.000 (-0.003) 


0.0053 (0.0049) 


0.006 


0.338 


0.217 


0.83 (0.83) 


2.20 (2.29) 


0.25 (0.26) 


0.003 ( 0.000) 


0.0008 (0.0009) 


0.006 


0.338 


0.340 


0.83 (0.83) 


2.16 (2.20) 


0.24 (0.25) 


-0.011 (-0.013) 


0.0008 (0.0009) 


0.006 


0.338 


0.507 


0.86 (0.87) 


2.03 (2.06) 


0.22 (0.22) 


-0.019 (-0.020) 


0.0012 (0.0014) 


0.007 


0.338 


0.743 


0.92 (0.92) 


1.87 (1.87) 


0.19 (0.19) 


-0.008 (-0.011) 


0.0030 (0.0034) 


0.010 



z 


Pt 
(GeV/c)^ 


(y) 


(GeV/c)^ 


{x) 




stat. uncertainty 


sys. uncertainty 


0.375 


0.005 


0.63 (0.64) 


2.50 (2.53) 


0.37 (0.36) 


0.006 ( 0.002) 


0.0032 (0.0029) 


0.007 


0.375 


0.026 


0.67 (0.67) 


3.13 (3.11) 


0.43 (0.43) 


0.018 (-0.010) 


0.0260 (0.0113) 


0.004 


0.375 


0.217 


0.84 (0.83) 


2.10 (2.19) 


0.23 (0.25) 


0.003 ( 0.001) 


0.0012 (0.0012) 


0.006 


0.375 


0.340 


0.84 (0.84) 


2.13 (2.19) 


0.24 (0.24) 


-0.007 (-0.010) 


0.0010 (0.0011) 


0.005 


0.375 


0.507 


0.86 (0.86) 


2.04 (2.07) 


0.22 (0.22) 


-0.015 (-0.018) 


0.0013 (0.0015) 


0.005 


0.375 


0.743 


0.92 (0.92) 


1.84 (1.88) 


0.19 (0.19) 


-0.011 (-0.014) 


0.0031 (0.0035) 


0.007 


0.413 


0.005 


0.66 (0.68) 


2.34 (2.38) 


0.33 (0.32) 


-0.011 (-0.011) 


0.0051 (0.0044) 


0.006 


0.413 


0.026 


0.65 (0.65) 


2.37 (2.37) 


0.34 (0.34) 


-0.003 (-0.002) 


0.1008 (0.0280) 


0.007 


0.413 


0.217 


0.87 (0.86) 


1.95 (1.99) 


0.21 (0.22) 


0.001 (-0.003) 


0.0023 (0.0023) 


0.004 


0.413 


0.340 


0.86 (0.86) 


2.17 (2.24) 


0.23 (0.24) 


-0.005 (-0.007) 


0.0012 (0.0014) 


0.007 


0.413 


0.507 


0.86 (0.87) 


2.04 (2.08) 


0.22 (0.22) 


-0.011 (-0.013) 


0.0015 (0.0017) 


0.007 


0.413 


0.743 


0.92 (0.92) 


1.85 (1.87) 


0.19 (0.19) 


-0.016 (-0.017) 


0.0032 (0.0036) 


0.004 


0.454 


0.005 


0.66 (0.68) 


2.35 (2.37) 


0.33 (0.32) 


0.000 (-0.006) 


0.0050 (0.0045) 


0.007 


0.454 


0.026 


0.68 (0.68) 


2.61 (2.59) 


0.35 (0.35) 


-0.007 ( 0.004) 


0.0401 (0.0141) 


0.004 


0.454 


0.340 


0.87 (0.87) 


2.02 (2.09) 


0.22 (0.22) 


-0.010 (-0.012) 


0.0016 (0.0018) 


0.005 


0.454 


0.507 


0.89 (0.89) 


1.98 (2.02) 


0.21 (0.21) 


-0.013 (-0.016) 


0.0020 (0.0023) 


0.004 


0.454 


0.743 


0.92 (0.92) 


1.85 (1.87) 


0.19 (0.19) 


-0.018 (-0.024) 


0.0035 (0.0041) 


0.005 


0.496 


0.005 


0.72 (0.72) 


2.63 (2.53) 


0.34 (0.32) 


0.008 (-0.010) 


0.0080 (0.0068) 


0.004 


0.496 


0.340 


0.87 (0.86) 


1.96 (2.00) 


0.21 (0.22) 


-0.011 (-0.011) 


0.0024 (0.0026) 


0.007 


0.496 


0.507 


0.92 (0.92) 


1.86 (1.88) 


0.19 (0.19) 


-0.014 (-0.017) 


0.0026 (0.0030) 


0.004 


0.496 


0.743 


0.93 (0.92) 


1.67 (1.68) 


0.17 (0.17) 


-0.011 (-0.014) 


0.0051 (0.0061) 


0.005 


0.540 


0.005 


0.71 (0.72) 


2.19 (2.20) 


0.28 (0.28) 


0.014 ( 0.011) 


0.0117 (0.0090) 


0.005 


0.540 


0.507 


0.92 (0.92) 


1.88 (1.91) 


0.19 (0.19) 


-0.021 (-0.024) 


0.0029 (0.0034) 


0.004 


0.540 


0.743 


0.94 (0.94) 


1.49 (1.49) 


0.15 (0.15) 


-0.019 (-0.042) 


0.0102 (0.0120) 


0.004 


0.586 


0.507 


0.92 (0.92) 


1.74 (1.74) 


0.18 (0.18) 


-0.016 (-0.018) 


0.0065 (0.0072) 


0.005 
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TABLE X: The extracted data on 7^4/(^2 + eHi) averaged over x and with their statistical and systematic uncertainties. 
The values given in brackets are obtained using the fit method. 



z 


Pt 
{GcV/cf 


(v) 


{GeV/cf 


(x) 


'Hil{H2+eHi) 


stat. uncertainty 


sys. uncertainty 


0.068 


0.026 


0.84 (0.84) 


2.19 (2.26) 


0.24 (0.25) 


0.001 ( 0.004) 


0.0017 (0.0021) 


0.018 


0.068 


0.065 


0.85 (0.85) 


2.15 (2.19) 


0.24 (0.24) 


0.037 ( 0.028) 


0.0013 (0.0018) 


0.011 


0.093 


0.026 


0.80 (0.79) 


2.23 (2.34) 


0.26 (0.28) 


0.006 ( 0.010) 


0.0011 (0.0012) 


0.013 


0.093 


0.065 


0.80 (0.80) 


2.29 (2.32) 


0.27 (0.27) 


0.012 ( 0.009) 


0.0007 (0.0009) 


0.015 


0.093 


0.129 


0.85 (0.85) 


2.11 (2.13) 


0.23 (0.23) 


0.041 ( 0.032) 


0.0013 (0.0016) 


0.014 


0.093 


0.217 


0.92 (0.92) 


1.86 (1.87) 


0.19 (0.19) 


0.025 ( 0.027) 


0.0066 (0.0076) 


0.007 


0.119 


0.026 


0.79 (0.78) 


2.18 (2.27) 


0.26 (0.27) 


0.009 ( 0.014) 


0.0016 (0.0017) 


0.013 


0.119 


0.065 


0.77 (0.77) 


2.31 (2.35) 


0.28 (0.29) 


0.005 ( 0.004) 


0.0008 (0.0009) 


0.013 


0.119 


0.129 


0.81 (0.81) 


2.24 (2.26) 


0.26 (0.26) 


0.028 ( 0.019) 


0.0010 (0.0012) 


0.014 


0.119 


0.217 


0.85 (0.85) 


2.01 (2.02) 


0.22 (0.22) 


0.051 ( 0.040) 


0.0020 (0.0024) 


0.014 


0.119 


0.340 


0.92 (0.92) 


1.68 (1.68) 


0.17 (0.17) 


0.048 ( 0.051) 


0.0142 (0.0157) 


0.008 


0.147 


0.005 


0.61 (0.59) 


2.75 (2.62) 


0.42 (0.41) 


0.024 (-0.034) 


0.0279 (0.0166) 


0.031 


0.147 


0.026 


0.87 (0.88) 


1.85 (1.90) 


0.20 (0.20) 


0.023 ( 0.027) 


0.0064 (0.0051) 


0.016 


0.147 


0.065 


0.75 (0.75) 


2.29 (2.34) 


0.29 (0.29) 


0.009 ( 0.010) 


0.0009 (0.0011) 


0.014 


0.147 


0.129 


0.78 (0.78) 


2.25 (2.28) 


0.27 (0.27) 


0.008 ( 0.005) 


0.0010 (0.0012) 


0.012 


0.147 


0.217 


0.82 (0.84) 


2.14 (2.43) 


0.24 (0.27) 


0.029 ( 0.023) 


0.0015 (0.0015) 


0.014 


0.147 


0.340 


0.91 (0.91) 


1.97 (1.99) 


0.20 (0.20) 


0.068 ( 0.060) 


0.0048 (0.0054) 


0.020 


0.175 


0.005 


0.59 (0.61) 


2.53 (2.60) 


0.39 (0.39) 


-0.008 (-0.011) 


0.0126 (0.0100) 


0.024 


0.175 


0.065 


0.77 (0.76) 


2.26 (2.33) 


0.28 (0.29) 


0.003 ( 0.004) 


0.0011 (0.0013) 


0.013 


0.175 


0.129 


0.78 (0.78) 


2.32 (2.35) 


0.28 (0.28) 


-0.004 (-0.006) 


0.0010 (0.0012) 


0.015 


0.175 


0.217 


0.81 (0.80) 


2.21 (2.23) 


0.26 (0.26) 


0.010 ( 0.007) 


0.0014 (0.0016) 


0.014 


0.175 


0.340 


0.85 (0.85) 


2.02 (2.03) 


0.22 (0.22) 


0.060 ( 0.049) 


0.0027 (0.0030) 


0.014 


0.175 


0.507 


0.92 (0.92) 


1.85 (1.86) 


0.19 (0.19) 


0.068 ( 0.062) 


0.0105 (0.0115) 


0.030 


0.205 


0.005 


0.61 (0.63) 


2.67 (2.71) 


0.40 (0.40) 


0.017 ( 0.018) 


0.0083 (0.0072) 


0.022 


0.205 


0.065 


0.77 (0.76) 


2.17 (2.25) 


0.26 (0.28) 


0.010 ( 0.012) 


0.0015 (0.0017) 


0.011 


0.205 


0.129 


0.77 (0.77) 


2.32 (2.35) 


0.28 (0.29) 


-0.008 (-0.011) 


0.0011 (0.0013) 


0.013 


0.205 


0.217 


0.80 (0.80) 


2.26 (2.28) 


0.26 (0.27) 


-0.007 (-0.010) 


0.0014 (0.0017) 


0.015 


0.205 


0.340 


0.84 (0.84) 


2.15 (2.17) 


0.24 (0.24) 


0.034 ( 0.031) 


0.0024 (0.0027) 


0.015 


0.205 


0.507 


0.92 (0.92) 


1.86 (1.87) 


0.19 (0.19) 


0.099 ( 0.084) 


0.0079 (0.0088) 


0.025 


0.236 


0.005 


0.61 (0.62) 


2.60 (2.62) 


0.39 (0.39) 


0.009 (-0.001) 


0.0063 (0.0060) 


0.023 


0.236 


0.065 


0.84 (0.83) 


1.87 (1.88) 


0.21 (0.21) 


0.006 ( 0.010) 


0.0044 (0.0044) 


0.009 


0.236 


0.129 


0.77 (0.77) 


2.31 (2.35) 


0.28 (0.29) 


0.001 (-0.002) 


0.0012 (0.0015) 


0.011 


0.236 


0.217 


0.78 (0.78) 


2.26 (2.28) 


0.27 (0.27) 


-0.008 (-0.011) 


0.0014 (0.0017) 


0.014 


0.236 


0.340 


0.81 (0.81) 


2.15 (2.16) 


0.25 (0.25) 


0.005 ( 0.007) 


0.0021 (0.0024) 


0.014 


0.236 


0.507 


0.87 (0.87) 


1.93 (1.93) 


0.20 (0.21) 


0.084 ( 0.078) 


0.0053 (0.0060) 


0.020 


0.236 


0.7 l:i 


0.92 (0.92) 


l.()7 (l.()8) 


0.17 (0.17) 


0.070 ( 0.071) 


0.()22() (().02r)2) 


().():', 1 


0.269 


0.005 


0.63 (0.64) 


2.66 (2.67) 


0.39 (0.38) 


0.007 ( 0.002) 


0.0063 (0.0062) 


0.022 


0.269 


0.129 


0.77 (0.77) 


2.23 (2.29) 


0.27 (0.28) 


0.011 ( 0.008) 


0.0015 (0.0018) 


0.012 


0.269 


0.217 


0.78 (0.78) 


2.26 (2.29) 


0.27 (0.27) 


-0.009 (-0.012) 


0.0015 (0.0018) 


0.013 


0.269 


o.;mo 


0.81 (0.81) 


2.16 (2.17) 


0.25 (0.25) 


0.007 ( 0.007) 


0.0021 (0.0024) 


0.013 


0.269 


0..j()7 


0.81 (0.81) 


2.02 (2.():i) 


0.22 (0.22) 


0.017 ( 0.012) 


0.()():i9 (0.1)1)1.")) 


0.020 


0.269 


0.743 


0.92 (0.92) 


1.79 (1.80) 


0.18 (0.18) 


0.088 ( 0.072) 


0.0162 (0.0179) 


0.021 


0.303 


0.005 


0.63 (0.64) 


2.51 (2.52) 


0.36 (0.36) 


0.013 ( 0.011) 


0.0069 (0.0069) 


0.021 


0.303 


0.129 


0.79 (0.78) 


2.06 (2.11) 


0.25 (0.25) 


0.010 ( 0.008) 


0.0023 (0.0026) 


0.010 


0.303 


0.217 


0.79 (0.79) 


2.27 (2.30) 


0.27 (0.27) 


-0.004 (-0.007) 


0.0017 (0.0021) 


0.013 


0.303 


0.340 


0.80 (0.80) 


2.17 (2.19) 


0.25 (0.25) 


0.005 ( 0.006) 


0.0022 (0.0025) 


0.013 


0.303 


0.507 


0.84 (0.84) 


2.02 (2.03) 


0.22 (0.22) 


0.034 ( 0.032) 


0.0039 (0.0044) 


0.022 


0.303 


0.743 


0.92 (0.92) 


1.86 (1.86) 


0.19 (0.19) 


0.103 ( 0.094) 


0.0127 (0.0144) 


0.028 


0.338 


0.005 


0.64 (0.65) 


2.50 (2.51) 


0.36 (0.36) 


-0.001 (-0.001) 


0.0067 (0.0069) 


0.021 


0.338 


0.026 


0.62 (0.62) 


2.66 (2.64) 


0.39 (0.39) 


-0.026 (-0.032) 


0.0351 (0.0203) 


0.022 


0.338 


0.129 


0.92 (0.92) 


1.74 (1.74) 


0.18 (0.18) 


0.001 (-0.004) 


0.0120 (0.0135) 


0.007 


0.338 


0.217 


0.79 (0.79) 


2.26 (2.31) 


0.27 (0.27) 


0.001 (-0.001) 


0.0020 (0.0024) 


0.010 


0.338 


0.340 


0.80 (0.80) 


2.17 (2.19) 


0.25 (0.26) 


0.005 ( 0.002) 


0.0023 (0.0028) 


0.010 


0.338 


0.507 


0.84 (0.84) 


2.02 (2.04) 


0.22 (0.23) 


0.020 ( 0.017) 


0.0039 (0.0046) 


0.010 


0.338 


0.743 


0.92 (0.92) 


1.86 (1.86) 


0.19 (0.19) 


0.087 ( 0.088) 


0.0124 (0.0138) 


0.015 



42 



z 


Pt 
{GeV/cf 


{y) 


(GeV/c)^ 


{x) 


ni/{n2 + eHi) 


stat. uncertainty 


sys. uncertainty 


0.375 


0.005 


0.64 (0.65) 


2.49 (2.49) 


0.36 (0.36) 


-0.007 (-0.007) 


0.0069 (0.0073) 


0.020 


0.375 


0.026 


0.67 (0.67) 


3.11 (3.11) 


0.43 (0.43) 


-0.004 (-0.036) 


0.0271 (0.0198) 


0.018 


0.375 


0.217 


0.80 (0.80) 


2.14 (2.20) 


0.25 (0.26) 


-0.004 (-0.005) 


0.0026 (0.0031) 


0.010 


0.375 


0.340 


0.81 (0.81) 


2.14 (2.17) 


0.25 (0.25) 


0.000 (-0.006) 


0.0027 (0.0033) 


0.009 


0.375 


0.507 


0.84 (0.84) 


2.04 (2.05) 


0.23 (0.23) 


0.013 ( 0.007) 


0.0042 (0.0049) 


0.009 


0.375 


0.743 


0.92 (0.92) 


1.86 (1.87) 


0.19 (0.19) 


0.067 ( 0.066) 


0.0124 (0.0143) 


0.010 


0.413 


0.005 


0.67 (0.68) 


2.34 (2.37) 


0.32 (0.32) 


-0.018 (-0.024) 


0.0105 (0.0106) 


0.015 


0.413 


0.026 


0.65 (0.65) 


2.37 (2.37) 


0.34 (0.34) 


-0.020 (-0.018) 


0.0851 (0.0357) 


0.018 


0.413 


0.217 


0.84 (0.83) 


1.97 (1.99) 


0.22 (0.22) 


0.002 (-0.006) 


0.0048 (0.0054) 


0.013 


0.413 


0.340 


0.84 (0.83) 


2.20 (2.24) 


0.24 (0.25) 


0.005 (-0.001) 


0.0033 (0.0040) 


0.010 


0.413 


0.507 


0.84 (0.84) 


2.03 (2.05) 


0.22 (0.23) 


0.003 (-0.003) 


0.0046 (0.0054) 


0.016 


0.413 


0.743 


0.92 (0.92) 


1.86 (1.88) 


0.19 (0.19) 


0.037 ( 0.036) 


0.0132 (0.0153) 


0.014 


0.454 


0.005 


0.68 (0.68) 


2.36 (2.36) 


0.32 (0.32) 


0.024 ( 0.011) 


0.0104 (0.0108) 


0.014 


0.454 


0.026 


0.67 (0.68) 


2.56 (2.58) 


0.35 (0.35) 


-0.053 (-0.031) 


0.0392 (0.0226) 


0.015 


0.454 


0.340 


0.84 (0.84) 


2.03 (2.06) 


0.23 (0.23) 


0.000 (-0.002) 


0.0041 (0.0050) 


0.017 


0.454 


0.507 


0.87 (0.87) 


1.95 (1.97) 


0.21 (0.21) 


0.005 (-0.002) 


0.0064 (0.0078) 


0.009 


0.454 


0.743 


0.92 (0.92) 


1.86 (1.87) 


0.19 (0.19) 


0.036 ( 0.016) 


0.0143 (0.0168) 


0.019 


0.496 


0.005 


0.72 (0.72) 


2.50 (2.49) 


0.32 (0.32) 


-0.017 (-0.037) 


0.0160 (0.0165) 


0.012 


0.496 


0.340 


0.83 (0.83) 


1.98 (2.01) 


0.22 (0.22) 


0.006 ( 0.009) 


0.0055 (0.0066) 


0.018 


0.496 


0.507 


0.90 (0.91) 


1.83 (1.85) 


0.19 (0.19) 


-0.008 (-0.007) 


0.0092 (0.0116) 


0.010 


0.496 


0.743 


0.92 (0.92) 


1.68 (1.68) 


0.17 (0.17) 


-0.014 (-0.005) 


0.0201 (0.0245) 


0.009 


0.540 


0.005 


0.71 (0.72) 


2.18 (2.20) 


0.28 (0.28) 


0.009 ( 0.005) 


0.0201 (0.0209) 


0.010 


0.540 


0.507 


0.92 (0.92) 


1.89 (1.90) 


0.19 (0.19) 


-0.006 (-0.005) 


0.0108 (0.0137) 


0.016 


0.540 


0.743 


0.94 (0.94) 


1.49 (1.49) 


0.15 (0.15) 


-0.006 (-0.081) 


0.0444 (0.0569) 


0.022 


0.586 


0.507 


0.92 (0.92) 


1.74 (1.74) 


0.18 (0.18) 


-0.013 (-0.010) 


0.0195 (0.0250) 


0.016 
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